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Espace vectoriel 


Introduction 


En géométrie plane, il est d'usage d’introduire la notion de vecteur. L'ensemble de ces 
vecteurs du plan est muni de deux opérations de base : l’addition des vecteurs, et la 
multiplication d’un vecteur par un nombre réel. Bien d’autres ensembles sont munis de 
deux opérations possédant des propriétés analogues à celles-ci (des ensembles de fonctions, 
de polynômes, etc). C’est cette structure commune, appelée structure d’espace vectoriel, 
que nous allons étudier. Compte tenu de l’origine concrète que nous avons indiquée, il ne 
faudra pas s’étonner du langage géométrique employé. De plus il est conseillé de s’appuyer 
dans une certaine mesure sur cette analogie pour nourrir son intuition concernant les 


structures introduites. 


1.1 La structure d’espace vectoriel 


Définition 1.1.1 Ün espace vectoriel E sur un corps commutatif K° (ap- 
pelé corps des scalaires) est un ensemble non vide muni d’une opération in- 
terne (addition qu'on notera +) telle que 


EV-1) E est un groupe abélien pour l’addition. 


et d’une opération externe, application de K x E dans E, qui à un couple 
(À, x) associe l’élément À.x (qu’on notera aussi quand aucune confusion n’en 
résulte Àx) telle que 


EV-2) À(x + y) = Àx + Ày. 
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EV-3) X{ux) = (Au)x. 
EV-4) (1 + ujx = Àx + ur. 
EV-5) 1x = x. 


Ex-1 Montrer que pour tout x € E et tout À € K, 0x = 0 et À0 = 0 
(attention par abus de notation il est d’usage de noter 0 le vecteur nul et 
aussi le zéro du corps de base; ce sont bien sûr des objets différents. Le 
lecteur prendra garde dans un premier temps à s’assurer qu’il est capable de 
bien distinguer ces deux objets dans les diverses formules écrites). 


Attention l'intuition géométrique à aussi ses limites, notamment quand le 
corps de base est de caractéristique non nulle. 


Ex-2 Montrer que si K est un corps de caractéristique 2, pour tout x € E, 
x + x = 0. Que se passe-t-il pour un corps de caractéristique p ? 


1.2 Premiers exemples d’espaces vectoriels 
On laisse le soin au lecteur de préciser les opérations et de vérifier les pro- 
priétés requises dans les exemples suivants. 


Ex-3 (Géométrie) L'ensemble des vecteurs de la géométrie plane ainsi que 
l’ensemble des vecteurs de la géométrie dans l’espace sont des espaces vecto- 
riels sur R. 


Ex-4 (Nombres complexes) C est un espace vectoriel sur R. € est un 
espace vectoriel sur C. 


Ex-5 (Nombres réels) R est un espace vectoriel sur Q. 


Ex-6 (Fonctions) L'ensemble des fonctions définies sur un intervalle 7 de 
R à valeurs dans R est un espace vectoriel sur R. 


Ex-7 L'ensemble des fonctions polynomiales de R dans R est un espace 
vectoriel sur R. L'ensemble des fonctions polynomiales de degré < n de R 
dans R est un espace vectoriel sur R. 


Ex-8 Tout corps commutatif K est un espace vectoriel sur lui même. 
Ex-9 K” est un espace vectoriel sur K. 


Ex-10 (Polynômes) Les polynômes à r variables et à coefficients dans un 
corps K forment un espace vectoriel sur K. 


Espaces vectoriels 9 


Ex-11 Les polynômes à r variables, à coefficients dans un corps K et de 
degré total < n, forment un espace vectoriel sur K'. 


Ex-12 Les polynômes à r variables, à coefficients dans un corps Æ et dont 
les degrés partiels par rapport à certaines des variables sont bornés, est un 
espace vectoriel sur K. 


Ex-13 Si Æ est un sous corps d’un corps commutatif X;, alors X, est un 
espace vectoriel sur K. 


Ex-14 Soit À un ensemble non vide, et F(A, E) l’ensemble des fonctions de 
À dans un espace vectoriel E sur le corps K. Alors F(A,E) est un espace 
vectoriel sur K. 


Ex-15 (Espace produit) Si E et F sont deux espaces vectoriels sur K, le 
produit E x Fest un espace vectoriel sur K. 


Ex-16 (Suites) Montrer que les suites d'éléments de À forment un espace 
vectoriel sur K. 


Ex-17 Que dire des suites réelles ? Des suites réelles qui convergent vers 0 ? 
Des suites réelles qui divergent vers +co ? Des suites réelles croissantes ? Des 
suites réelles convergentes ? 


Ex-18 (Algèbre de Boole) Soit À un ensemble non vide et P(A) l’ensemble 
des parties de A. On définit sur P(A) la somme de deux parties comme étant 
leur différence symétrique (U + V = UAV = {(U-—-V)U(V —U)). De plus 
on définit le produit d’un élément de P(A) par un élément du corps à deux 
éléments {0,1} par O.U = Ÿ et 1.U = U. Montrer que dans ces conditions 
P(A) est un espace vectoriel sur le corps à deux éléments. 

Remarque : Si de plus on tient compte de l’opération multiplication définie 
comme étant l’intersection de deux parties, on a une structure d’algèbre de 
Boole . 
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Sous espaces vectoriels 


Introduction 


Comme toujours, lorsqu'on définit une structure, on introduit la notion de sous structure. 
Ainsi nous étudions maintenant la notion de sous espace vectoriel. Beaucoup d’ensembles 
intéressants sont des sous espaces vectoriels d'espaces plus gros. Dans ce cas pour montrer 
que ce sont des espaces vectoriels il n’est pas nécessaire de vérifier toutes les propriétés de 
l’addition et de la multiplication externe. La stabilité par rapport à ces opérations suffit. 


2.1 La sous structure de sous espace vectoriel 


Définition 2.1.1 Soit E un espace vectoriel sur un corps K. Un sous en- 
semble F de E est un sous espace vectoriel de E si c’est un espace vectoriel 
pour les deux opérations induites par E sur F. 


Ex-1 Montrer qu’un sous ensemble non vide F d’un espace vectoriel Æ sur 
un corps commutatif X est un sous espace vectoriel de Æ si et seulement si 
Va,yeF, x+yer, 

VAeK,VreF, rer. 

(on dit alors que F est stable par rapport à ces opérations). 
Ceci donne un moyen commode de montrer que des ensembles sont munis 


d’une structure d’espace vectoriel. 


LL 
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Ex-2 Montrer que l'intersection d’une famille non vide de sous espaces 
vectoriels d’un espace Æ est un sous espace de E,. 


Ex-3 Peut on dire que la réunion d’une famille de sous espaces vectoriels 
de E est un sous espace vectoriel de E ? 


Définition 2.1.2 Soit À un sous ensemble de l’espace vectoriel E. Considérons 
la famille de tous les sous espaces de E qui contiennent À. Cette famille est 
non vide puisque E en fait partie. L’intersection des sous espaces de cette 
famille est appelée le sous espace engendré par A. On note [A] ce sous 
espace. 


Ex-4 Le sous espace engendré par À est le plus petit sous espace de E (au 
sens de l'inclusion des ensembles) contenant A. 


2.2 Exemples de sous espaces vectoriels 


Ex-5 Montrer que l’ensemble C'(R) des fonctions continues sur R à valeurs 
dans R est un sous espace vectoriel de l’espace de toutes les fonctions définies 
sur R à valeurs dans R. 


Ex-6 Montrer que l’ensemble des fonctions dérivables sur R à valeurs dans 
R est un sous espace vectoriel de l’espace de toutes les fonctions définies sur 
R à valeurs dans R. Que peut-on en dire par rapport à C(R) ? 


Ex-7 Montrer que l’ensemble F,(A,E) des fonctions définies sur un en- 
semble non vide À, à valeurs dans un espace vectoriel E, et qui s’annulent 
en un point fixé a € À, est un sous espace de l’espace F(A, E) des fonctions 
de À dans E. 


Ex-8 Montrer que l’ensemble des fonctions P(R) de R dans R, qui sont 
paires est un sous espace vectoriel de l’espace F(R, R) de toutes les fonctions 
de R dans R. Est-ce aussi vrai pour les fonctions impaires ? 


Ex-9 Soit [a, b] un intervalle fixé de R. Une fonction f de [a, b] dans R est 
dite en escalier s’il existe un partage a = ï0 < æ1 < +-: < xn = b de [a, b] tel 
que f soit constante sur chaque intervalle [x,x;,1[. Montrer que l’ensemble 
des fonctions en escalier est un sous espace vectoriel de l’espace F([a, b], R) 
de toutes les fonctions de {a, b] dans R. 
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Ex-10 Montrer que les suites de nombres complexes qui convergent vers 0 
forment un sous espace vectoriel de l’espace de toutes les suites de nombres 
complexes. 


Ex-11 Soit F l’espace des fonctions continues sur un intervalle [a, b], affines 
(i.e. de la forme ax+b) par morceaux (ï.e. il existe un partage a = xo < t1 < 

+ < Zn = b de [a,b] tel que f soit affine sur chaque intervalle [x;x;,:]). 
Montrer que F est un sous espace de F({[a, b|, R). 
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Homomorphismes d’espaces 
vectoriels 


Introduction 


Nous étudions maintenant les applications d’un espace vectoriel dans un autre (sur un 
même corps) qui conservent la structure. Ces applications, dites applications linéaires, 
sont centrales dans la théorie. Nous reviendrons dans la plupart des chapitres qui suivent 
à leur étude approfondie. 


3.1 Application linéaire - Noyau - Image 


Définition 3.1.1 E et F étant deux espace vectoriels sur le même corps K, 
un homomorphisme de E dans F (ou encore une application linéaire 
de E dans F) est une application f de E dans F telle que 


Hom-1) f(x + y) = f(x) + f(y). 

Hom-2) f(x) = Àf(x). 
Dans le cas où F = E on dit que f est un endomorphisme de E, ou encore 
un opérateur de E. 
Ex-1 Montrer que si f est une application linéaire f(0) = 0. 


Ex-2 Montrer que l’ensemble des éléments de Æ dont l’image est 0 est un 
sous espace vectoriel de Æ. 
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Ex-2 Montrer que si f est injective alors 0 est le seul élément de Æ dont 
l’image soit 0. 


Ex-3 Réciproquement, montrer que si 0 est le seul élément de E dont l’image 
soit 0 alors f est injective. 


Ex-4 Montrer que l’image de f (c’est-à-dire f(Æ)), est un sous espace vec- 
toriel de F. On note Zm(f) cette image. 


Définition 3.1.2 Si f est un homomorphisme de E dans F le noyau de f, 
noté Ker(f), est le sous espace vectoriel de E constitué des éléments de E 
qui ont pour image 0. 


Le noyau de f, (Ker(f)), et l’image de f, (1m(f)), sont deux objets très 
importants pour l’étude d’un homomorphisme f. Compte tenu de l’exercice 
2, on peut énoncer le résultat 


Théorème 3.1.1 Un homomorphisme f est injectif si et seulement si son 
noyau est {0}. Dans ce cas f est un isomorphisme de E sur son image notée 
Im(f) ou encore f(E). 


3.2 Premières propriétés des applications linéaires 


Dans la suite K désigne un corps commutatif. 
Ex-5 Citer des applications linéaires et des applications non linéaires. 


Ex-6 Montrer que dans l’espace vectoriel des vecteurs de la géométrie dans 
l’espace, la projection sur un plan parallèlement à une droite et la projection 
sur une droite parallèlement à un plan sont des applications linéaires dont 
on determinera le noyau et l’image. 

Soient 0 < n < m deux entiers. Soit f l’application de K”* dans K7 qui à 
(t1,---, Zm) fait correspondre (x1,---,x,). Montrer que f est un homomor- 
phisme (application linéaire) dont on déterminera le noyau et l’image. 


Ex-7 Soit C!(R) l’espace vectoriel des fonctions dérivables sur R et ayant 
une dérivée continue (fonctions de classe C1) et C(R) l’espace des fonctions 
continues sur R. Montrer que la dérivation D est un homomorphisme de 
CT(R) dans C(R). Quels sont son noyau et son image ? 
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Ex-8 Soit K[X] l’espace vectoriel des polynômes à une variable à coeffi- 
cients dans K. Soit T l'application de K[X] dans lui même qui à un po- 
lynôme P(X) fait correspondre le polynôme P(X + 1). Montrer que T est 


un endomorphisme. Quels sont le noyau de T et son image ? 


Ex-9 Soit K,[X] l’espace vectoriel des polynômes à une variable à coeffi- 
cients dans X et de degé < n. Soit D l'application de K,[X] dans lui même 
qui à un polynôme associe son polynôme dérivé. Montrer que D est un en- 
domorphisme. Quels sont son noyau et son image ? 


Ex-10 Soit X[X1] l’espace vectoriel des polynômes à une variable à coeffi- 
cients dans X. Soit F(K) l’espace des fonctions de K dans lui même. Montrer 
que l’application 6 qui à un polynôme P(X) fait correspondre la fonction po- 
lynômiale x + P(x) est un homomorphisme. On suppose le corps K infini, 
quel est le noyau de 4? 

Supposons maintenant que Æ soit le corps à deux éléments, quel est le noyau 
de ®? 


Ex-11 Soit K[X] l’espace vectoriel des polynômes à une variable à coeffi- 
cients dans K et soit D(X) un polynôme non nul. On considère l’application 
de XT[X] dans lui même qui à tout polynôme P(X) associe le reste de la 
division euclidienne de P(X) par D(X). Montrer que cette application est 
un endomorphisme. Quels sont son noyau et son image ? 


Ex-12 Soit X,[X] l’espace vectoriel des polynômes à une variable à coef- 
ficients dans K et de degé < n. On fixe un point a € K et on considère 
l'application de K,[X] dans K qui à tout polynôme P(X) associe sa valeur 
P(a). Montrer que cette application est un homomorphisme. Quels sont son 
noyau et son image ? 


Ex-13 Soit C(R) l’espace des fonctions continues sur R à valeurs dans R. 
Fixons deux points a et b dans R. Soit  l’application de C'(R) dans R définie 
par 


b 
1 = | ot 
a 
Montrer que cette application est un homomorphisme. Quelle est son image ? 


Ex-14 Fixons a dans R. Soit f l'application de R dans R° qui à x associe 
(x, ax). Montrer que f est une application linéaire. Quels sont son noyau et 
son image ? 
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Ex-15 Fixons a, b, c, d dans R. Soit f l’application de R° dans R? qui à tout 
(x, y) associe (ax + by, cx + dy). Montrer que f est un endomorphisme. Quels 
sont son noyau et son image ? 


Le cas où f est une application linéaire de Æ dans le corps de base K est 
particulier. On dit que f est une forme linéaire. Ce cas, par son importance, 
motive une étude spécifique ultérieure. 


3.3 Structure de l’ensemble des applications 
linéaires 


Ex-16 Soient ÆE et F deux espaces vectoriels sur le même corps X. Notons 
£(E, F) l’ensemble des homomorphismes de Æ dans F. Si f et g sont des 
éléments de £(E, F) et si À est un élément de X, on définit f +g, Àf comme 
on le fait habituellement pour les fonctions, c’est-à-dire par 


(f + g)(x) = f(x) + g(x), 
(AF)(x) = A f(x). 
Montrer que muni de ces opérations, £L(E, F) est un espace vectoriel sur X. 


Ex-17 Soit E,F,G trois espaces vectoriels sur le même corps K. Si f est 
une application linéaire de Æ dans F et g une application linéaire de F dans 
G, montrer que go f est une application linéaire de E dans G. 


Ex-18 On considère maintenant l’espace £L(E) des endomorphismes (des 
opérateurs) de E. En plus de l’addition et de la multiplication par un sca- 
laire définies dans un exercice précédent, on considère aussi l’opération de 
composition des applications f o g. Montrer que ces opérations confèrent à 
£(E) une structure d’algèbre (non commutative). 

Montrer que si f et g éléments de L(E) vérifient f og = I (I est l'application 
identité), alors f et g ont des inverses. 
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Combinaisons linéaires 


Introduction 


Comme dans le cas du plan géométrique ou de l’espace, il est important de savoir dans 
quelles conditions on peut décomposer un vecteur sur d’autres vecteurs et dans quels cas 
cette décomposition est unique. Ces questions conduisent aux notions de familles libres, 
de familles génératrices et enfin à la notion de base. En examinant les exemples donnés, 
le lecteur pourra faire un inventaire (non exhaustif) des méthodes utilisées pour montrer 
que des familles sont libres, que des familles sont génératrices, que des familles sont des 
bases. 


4.1 Partie engendrée par une famille de vec- 
teurs 
On rappelle qu’une famille non vide d’éléments d’un ensemble À est une 


application 4 d’un ensemble non vide 7 (appelé ensemble des indices) dans 
A. La famille est alors notée (x;)ier où x; = O(i) (€ À). 


Définition 4.1.1 Soit E un espace vectoriel sur un corps K. Soit S = (xi)ier 
une famille non vide de E. On appelle combinaison linéaire d'éléments 
de S toute somme de la forme 


iel 


où seuls un nombre fini de coefficients À; sont non nuls. 
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Ex-1 Montrer que l’ensemble des combinaisons linéaires d'éléments de S 
est un sous espace vectoriel de Æ, et que ce sous espace est le sous espace 
engendré par S. On notera [S] ce sous espace. 


Définition 4.1.2 La famille S d'éléments de E est dite génératrice si 
[S] = Æ. {Autrement dit tout élément de E s'exprime comme combinaison 
linéaire des éléments de S). 


Définition 4.1.3 La famille S d'éléments de E est dite libre s4 tout élément 
de [S] s'écrit de manière unique comme combinaison linéaire des éléments 
de S (Les éléments de S sont alors dits linéairement indépendants). Si 
la famille S n’est pas libre, on dit qu’elle est liée {les éléments de S sont 
alors linéairement dépendants). 


Remarque : attention, on travaille ici sur des familles, et non simplement 
sur des parties. Par exemple (x:1,41,x%2) est une famille de trois éléments ; si 
T1 £ Lo, la partie qui lui est associée est {x1,x2}. 


Ex-2 Montrer que S'est une famille libre si et seulement si on a l’implication 
suivante 
ie 

Ex-3 Montrer que la famille S'est liée, si et seulement s’il existe un vecteur 
de $ qui est combinaison linéaire des autres vecteurs de S$. 
Montrer que si la famille S contient 0, alors elle est liée. 
La réciproque est elle vraie ? 
Montrer que si la famille S contient deux vecteurs identiques, alors elle est 
liée. 
Ex-4 Soit Æ un espace vectoriel non réduit à {0}. Montrer l’équivalence 
entre les trois propriétés suivantes 

1) S'est une famille libre maximale au sens de l’inclusion dans l’ensemble 
des parties libres de E. 

2) S est une famille génératrice minimale au sens de l’inclusion dans 
l’ensemble des parties génératrices de E. 

3) S est une famille à la fois libre et génératrice. 
Que se passe-t-il dans le cas particulier très simple d’un espace vectoriel 
réduit à {0}? 
Définition 4.1.4 Une famille libre et génératrice de E est appelée une base 
de E. 
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4.2 Exemples de familles de vecteurs 


Ex-5 Donner des exemples de la géométrie plane et de la géométrie dans 
l’espace de familles libres et non génératrices, de familles génératrices et non 
libres, de familles libres et génératrices. 


Ex-6 Soit À un ensemble non vide, E un espace vectoriel. On note F(A, E) 
l’espace des fonctions de À dans E. Soit S la partie de F(A, FE) constituée des 
fonctions constantes et des fonctions qui s’annulent en un point fixé 0 € À. 
Déterminer [S]. 


Ex-7 Soit [a, b] un intervalle de R. Notons F([a, b], R) l’espace des fonctions 
de [a, b] dans R. On considére le sous ensemble S de F{({a, b], R) défini par 
S={xmx-u+ir-ulue la,bl}. 


Déterminer [S]. 
Déterminer [S U {1}]. 


Ex-8 Soit (;);>1 une suite de nombres réels dont les termes sont tous 
distincts. Montrer que la famille de fonctions 


(es) 


est libre. 
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Bases - Dimension 


Introduction 


La notion de base est centrale pour l’étude des espaces vectoriels de dimension finie. Tout 
espace vectoriel admet une base, mais en dimension infinie, dans la plupart des cas il est 
impossible d’en construire une explicitement. Comme nous le verrons dans la suite, on est 
souvent amené à considérer plusieurs bases sur un même espace, et regarder alors quels 


sont les liens entre les composantes des vecteurs dans les diverses décompositions. 


5.1 Bases d’un espace vectoriel 


Définition 5.1.1 Une base d’un espace vectoriel E est une famille (e;)ier 
d'éléments de E qui est à la fois libre et génératrice. 


Définition 5.1.2 Un espace E est dit de dimension finie s’il possède une 
famille génératrice finie. 


Ex-1 Montrer que tout espace vectoriel de dimension finie non réduit à 
{0} admet une base ayant un nombre fini d'éléments. (En fait tout espace 
vectoriel non réduit à {0} admet une base, mais la démonstration quand on 
n’est pas en dimension finie n’est pas du niveau de ce cours). 


Le résultat qui suit (théorème de l’échange), est fondamental dans l’étude 
des bases des espaces vectoriels. 
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Ex-2 Soit E un espace vectoriel, F une famille libre de Æ et G une famille 
génératrice de Æ. Montrer que si [F] # E alors il existe x € G tel que FÜ{x} 
soit une famille libre, et que dans ces conditions [FT C [FU {x}]. 


Ex-3 (Théorème de la base incomplète). Soit Æ un espace vectoriel 
de dimension finie, F une famille libre de Æ et G une famille génératrice 
de E. Montrer qu’on peut compléter la famille F en une base de E en lui 
adjoignant des vecteurs de G. 


Ex-4 Montrer que dans un espace vectoriel de dimension finie, toutes les 
bases ont le même nombre d’éléments. 


Définition 5.1.3 Si E est un espace de dimension finie nous appellerons 
dimension de E et noterons dim(Æ) le nombre d'éléments d’une base de 
E (toutes les bases ont le même nombre d'éléments) si E Æ {0}, et O si 
PEAUX 


Remarque : en dimension finie n, on indexe en général les éléments d’une 
base par {1,2,---,n} (ou par {0,1,---,n—1}). 


Ex-5 Soit E un espace vectoriel de dimension finie n. Montrer que (e;)1<i<n 
est une base de Æ si et seulement si e; Æ 0 et pour tout 2<i<ne; & 
ex, ire CEE 


Définition 5.1.4 Soit E un espace vectoriel, (s;);er une famille de vecteurs 
de E telle que le sous espace S = [silier engendré par les vecteurs s; soit de 
dimension finie. On appelle rang de la famille (5;);er, la dimension de S. 


Remarques : 


1) Sur les méthodes utilisées en dimension finie. En dimension fi- 
nie, si on connaît la dimension n de l’espace, pour montrer qu’une famille 
de n éléments est une base, il suffit de montrer qu’elle est soit libre, soit 
génératrice. Remarquons qu’un espace de dimension n sur un corps fini à q 
éléments a g” éléments ; cette remarque permet de trouver la dimension n de 
l’espace si on connaît son cardinal. 


2) Suites et fonctions. Une fonction f d’un ensemble fini à n éléments 
dans un corps K,, outre sa représentation fonctionnelle, peut être définie aussi 
par ses valeurs, c’est-à-dire par un n-uple (a,---,a,). Ce n-uple peut être 
aussi considéré comme la suite des coefficients d’un polynôme a + a2X + 
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+ anX"" 1, Ces trois représentations d’un même objet peuvent être utiles. 
De plus se donner une base d’un espace de dimension n, revient à représenter 
tout élément de cet espace par ses composantes, c’est-à-dire par un n-uple, 
ce qui donne aux trois représentations précédentes une grande importance. 


3) Droites et hyperplans. Dans un espace vectoriel Æ de dimension n, 
un sous espace vectoriel de dimension 1 est appelé droïte vectorielle (ou 
plus brièvement droiïte), un sous espace vectoriel de dimension n — 1 est 
appelé hyperplan vectoriel (ou plus brièvement hyperplan). 


5.2 Exemples de bases 
Ex-6 Montrer que dans K”, la famille (e;)1<i<, où 
Ei — (0,---,0,1,0,---,0) 


(le 1 étant la i* composante) est une base. 


Ex-7 Montrer que la famille (1, X, X?,---) est une base de X[X] et que la 
famille (1,X,X2,---, X”) est une base de l’espace K,[X] des polynômes de 
degré < n. 


Ex-8 Soit F(A, K) l’espace des fonctions définies sur un ensemble fini À à 
valeurs dans le corps K. Pour tout a € À, notons e, la fonction définie par 
e(a) = 1 et e(x) = 0 si x £ a. Montrer que la famille (e,),e4 est une base. 
Quelles sont les composantes d’une fonction f sur cette base ? 


Ex-9 Dans l’espace K,[X] des polynômes de degré < n on considère pour 
tout 0 < i < n un polynôme P;(X) de degré exactement à. Montrer que la 
famille (P(X)),<<, est une base de K,[XT. 


Ex-10 Soit E un espace vectoriel de dimension n et (e;)1<;<, une base de 
E. Pour tout 1 < j < n on pose 


J = ÿ Qi, j€i 


1<i<ÿ 


avec à; Z 0. Montrer que la famille (f;)1<;<n est une base de E. 


Ex-11 Soit Æ un espace vectoriel de dimension n sur le corps fini F, à 
q éléments. Quel est le nombre d’éléments de E ? Quel est le nombre de 
familles libres de Æ ayant k éléments ? Quel est le nombre de bases de E ? 
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Combien y a-t-il de sous espaces de dimension 4 dans E ? Quel est le nombre 
de directions de droites dans E ? 


Ex-12 Soit B, = F({0,1}",{0,1}) l’espace (sur le corps à deux éléments) 
des fonctions booléennes de n variables booléennes. Quel est le nombre d'éléments 
de B, ? En déduire la dimension de B,. Pour tout a = (&,--:,an) € {0,1}" 

on définit la fonction 


ealtis tn) = El (1 + x + &). 


1<i<n 


Montrer que la famille (ec)acgo1}" est une base de B, (Comparer avec l’exer- 
cice 8). 

Quelles sont les composantes d’une fonction f sur cette base ? 

Montrer que toutes les fonctions de B, sont des fonctions polynomiales en n 
variables de degré total au plus n et de degré partiel par rapport à chaque 
variable au plus 1. 

Trouver une base de B,, constituée de fonctions monomiales. 

Traiter en détail l'exemple n = 3. 


Ex-13 Soient E et F deux espaces vectoriels de dimension finie. Construire 
une base de l’espace L(E, F) des applications linéaires de E dans F. Quel 
est la dimension de cet espace ? 


Chapitre 6 


Somme directe de sous espaces 


Introduction 


La notion de somme de sous espaces, permet une décomposition plus grossière des vecteurs 
que la notion de base. On obtient par exemple des décompositions du type : toute fonction 


réelle est de manière unique somme d’une fonction paire et d’une fonction impaire. 


6.1 Décomposition d’espaces vectoriels 


Définition 6.1.1 Soient EF; et E2 deux sous espaces vectoriels d’un espace 
E. La somme de ces deux sous espaces, notée FE; + E:, est l’ensemble des 
T1 + To OÙ t1 € En et 2 € E. 


Ex-1 Montrer que la somme ÆE; + F2 est un sous espace vectoriel de E. 


Ex-2 Soit F = EF; + F2. Montrer que pour que tout élément de F s’écrive 
de manière unique sous la forme x, + z2 avec 1 € FE et 2 € Eo, il faut et 
il suffit que E1 N E2 = {0}. 


Définition 6.1.2 Si E; et F2 vérifient E1 N E2 = {0}, on dit que la somme 
FE, + FE, est une somme directe et on la note F1 @ Es. 


Ex-3 Soient Æ: et F2 deux sous espaces de dimension finie d’un espace E 
tels que Æ1 N E2 = {0}. Soit (u:,--:,u,) une base de Æï et (w1,---,v,) une 


base de Æ2. Montrer que (u1,---,u,,01,::+,0) est une base de Æ1 @ E2. 


2L 
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On en conclut que dim(Æ1 @ E2) = dim(E1) + dim(ÆE>). 


Ex-4 Soit E un espace vectoriel et (e;);er une base de Æ. Considérons une 
partition 7? = JU K de I. Montrer que si on pose Æ1 = [eilies et Eo = [eiliex, 
alors E = FE; e FE. 

Ex-5 Soit Æ un espace de dimension finie et E; un sous espace de E. Montrer 


qu'il existe un sous espace F2 de E tel que E = FE; @ F3. Ce sous espace Æ2 
est-il unique ? 

Ex-6 Si E: et E2 sont deux sous espaces de dimension finie de Æ, montrer 
que dim(Æ1 + E2) = dim(ÆE:) + dim(E>) — dim(E: N E2). 


Définition 6.1.3 Soit E un espace de dimension finie et E1 un sous es- 
pace de E. Un sous espace E2 de E tel que E = EF; @ E2 est appelé un 
supplémentaire de E; dans E. 


Ex-7 Montrer que tous les supplémentaires de E; ont la même dimension. 


Définition 6.1.4 En dimension finie, la dimension d’un supplémentaire de 
E\ est appelée la codimension de E1. et notée codim(E:). 


Ex-8 Soit E un espace de dimension finie dont on fixe une base (e1,---,e,). 
Soit C un sous espace vectoriel de dimension 4 de Æ. Montrer qu’il existe 
une sous famille de n — k vecteurs (e;,,---,e;,_x%) de la base (e1,---,e,) telle 
que 

E=Ce lee ses 2]. 
Pour tout j € {1,-:-,n} \{ü,:::,i,-,} on décompose e; en la somme d’un 
vecteur c; € C et d’un vecteur f; € [e;,,:-:,e;, ,] 

€; — Cj + re 

Montrer que la famille (c;) ainsi obtenue est une base de C. 
Comment sont les composantes des vecteurs e; où j € {1,---,n}\{i,"--,in-x} 
dans la base (gi)1<ien Où gi = ei si à € {ü,** nr} et gi = & sinon? 


Définition 6.1.5 Soit E un espace vectoriel, décomposé en somme directe 
de deux de ses sous espaces E; et EF: 


E=E, @ E:. 


La projection p sur E; parallèlement à E2 est l'application qui à tout x de E 
fait correspondre sa composante x, sur FE, dans la décomposition x = x1 +2 
où x1 € En et %2 € Fo. 
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Ex-9 Montrer que la projection p ainsi définie est une application linéaire 
vérifiant p? = p. Quel est le noyau de p ? Quelle est l’image de p? 


Ex-10 Soit E un espace vectoriel de dimension finie et f un endomorphisme 
de E tel que f? = f. Montrer que f est une projection. 


6.2 Exemples de décompositions 


Ex-11 Donner des exemples et des contre exemples dans le plan et dans 
l’espace (faire des dessins). 


Ex-12 Montrer que C =R@ fi]. 


Ex-13 Soient E et F deux espaces vectoriels sur le même corps. On considère 
les deux sous espaces ÆE1 et F, de E x F définis par Æ1 = E x {0} et 
F, = {0} x F. Montrer que E x F = FE @ F1. 


Ex-14 Soit F(R,R) l’espace de toutes les fonctions de R dans R. On note 
P le sous espace constitué des fonctions paires et Z le sous espace constitué 
des fonctions impaires. 

Montrer que F(R,R) = P @Z. 


Ex-15 Soit F(R, R) l’espace de toutes les fonctions de R dans R. On note 
F, le sous espace des fonctions nulles au point a, et C le sous espace des 
fonctions constantes. 

Montrer que F(R,R) = F, @C. 


x 


Ex-16 Soit K[X] l’espace des polynômes à une variable et à coefficients 
dans le corps K. Fixons N(X) € K[X] et notons n (supposé > 1) le degré 
de N(X). Soit E le sous espace de X[X] constitué des multiples A(X)N(X) 
de N(X). Montrer que 


KIX] = K,_1|X]@E, 


où Kh_1[X] est le sous espace constitué des polynômes de degré < n — 1. 


Ex-17 Soit Æ un espace vectoriel de dimension n, H un hyperplan vectoriel 
de Æ, D une droite vectorielle qui n’est pas contenue dans H. Montrer que 
E=HeD. 
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6.3 Somme de plusieurs sous espaces 


Nous étendons maintenant la notion de somme directe à un nombre fini 
quelconque de sous espaces. 


Ex-18 Montrer que la somme de deux sous espaces vectoriels est associative. 
Ceci nous permet donc d'envisager de généraliser la notion de somme à un 
nombre fini de sous espaces. 


Définition 6.3.1 Soit (Ei)1<;<, une famille finie de s sous espaces vectoriels 
d’un espace vectoriel E. La somme de ces s sous espaces notée É5 +---+E, 
est l’ensemble des x1 +---+x, où x; € E;. 


Ex-19 Montrer que la somme Æ, +---+ FE, est un sous espace vectoriel de 
FE. 


Ex-20 Montrer que pour que tout élément de E1 -:-+EÆ, s’écrive de manière 
unique sous la forme x, +-:-+ x, où x; € E,, il faut et il suffit que pour 
chaque 1 < à < s on ait 


FE; N (Ei +. + E;-_; + Er + cdots + E;) — {0}. 
Définition 6.3.2 Si la famille (Ei)1<;<s vérifie pour chaque 1 <i<s 
ES Er SE ie be PRE) 0} 


on dit que la somme El +--..+EÆ, est directe. On la note 


li<s 


Ex-21 Soit 
ESEÇG-.CE, 


une suite de sous espaces vectoriels de dimension finie d’un espace vectoriel 
E. Montrer qu’il existe une suite (F;)1<i<n de sous espaces vectoriels de E 


telle que pour tout 1<j<n 
1<i<ÿ 


Ex-22 Généraliser les notions de somme et de somme directe à une infinité 
de sous espaces. 


Chapitre 7 


Projections 


Introduction 


La notion de projection est une notion très importante liée à la décomposition en somme 
directe. On peut intuitivement se faire une idée de ce que représente une projection en se 
plaçant en géométrie à 3 dimensions sur R et en pensant par exemple à la projection sur 
un plan parallèlement à une droite. 


7.1 Projection et somme directe 


Définition 7.1.1 Soit E un espace vectoriel, décomposé en somme directe 
de deux de ses sous espaces E; et E2 


E=E, @ E:. 


La projection p sur EF; parallèlement à E2 est l’application de E dans 
E qui à tout x de E fait correspondre sa composante x: sur FE; dans la 
décomposition x = x1 + 2 où x1 € E1 et x2 € Eo. 


Ex-1 Montrer que la projection p ainsi définie est une application linéaire 
vérifiant p? = p. Quel est le noyau de p? Quelle est l’image de p. 

Ex-2 Soit E un espace vectoriel de dimension finie et f un endomorphisme 
de Æ tel que f? = f. Montrer que f est une projection. 
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Ex-3 Soit E un espace vectoriel de dimension finie, E; € F2 deux sous es- 
paces de Æ. Soit F> un supplémentaire de F2 dans E et F;, un supplémentaire 
de FE; dans F2. Montrer que la somme F; + F, est directe. 

Soit p2 la projection de Æ sur ÆE2 parallèlement à F2 et p. la projection de 
FE; sur FE; parallèlement à F,. Montrer que p = p1 © p2 est la projection de E 
sur ÆE; parallèlement à F, @ F2. 


Ex-4 Montrer que si deux projections p1 et p2 commutent alors p1 o p2 est 
une projection. 
Montrer que 

Im(p1 0 pa) = Im(p1) N Im(p2) 


et que 
Ker(p: o p2) = Ker(p1) + Ker(p2). 


Construire des exemples de projections qui commutent. 


Chapitre 8 


Formes linéaires 


Introduction 


Les formes linéaires sont des applications linéaires particulières : l’espace d’arrivée est de 
dimension 1, c’est le corps de base K. Ceci motive une étude spécifique. En effet, il est 
clair par exemple que l’étude d’une fonction de R” dans R” se ramène à l’étude de m 
fonctions de R° dans R (quand on veut étudier un signal où il y a n entrées et m sorties, 
on peut étudier chaque sortie en fonction des n entrées) ; en particulier on peut considérer 
que l’étude d’une application linéaire de R°* dans R”* peut passer par l’étude de m formes 
linéaires sur R”. Cette façon de voir, même si elle ne constitue pas le seul angle d’attaque, 
motive une étude particulière du cas des formes linéaires. De plus, du fait que m = 1, les 
résultats s'expriment souvent de manière plus agréable et plus instructive. 


8.1 Formes linéaires - Espace dual 


Définition 8.1.1 Soit E un espace vectoriel sur un corps K. Une forme 
linéaire sur E est une application linéaire de E dans K. 


Ex-1 Montrer que si f est une forme linéaire non nulle alors Im(f) = K 
(en particulier la dimension de l’image de f est 1). 


Définition 8.1.2 Le dual de E, que nous noterons E*, est l’ensemble des 
formes linéaires sur E. 
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Ex-2 Montrer que E*, muni de l’addition des formes linéaires et de la 
multiplication par un scalaire, est un espace vectoriel sur K. 


Ex-3 Soit f une forme linéaire non nulle sur un espace vectoriel Æ de 
dimension finie n. Montrer que le noyau de f est un hyperplan de E. 


Ex-4 Réciproquement, montrer que tout hyperplan de E est le noyau d’une 
forme linéaire. 


Ex-5 Soit E un espace vectoriel de dimension finie n muni d’une base 
(ei)icien de E. On considère dans £* la famille (e*)1<;<n définie par eÿ(ei) = 
d;,; (on rappelle que la notation de Kronecker 6; ; désigne 1 sii = j et 0 sinon). 
Que peut on dire de la dimension de E*. Montrer que la famille (ef)1<5<n est 
une base de Æ*. La dimension de E* est donc n. 


Définition 8.1.3 La base (eï)1<;j<n de E* est appelée la base duale de la 
base (ei)1<j<n- 


Définition 8.1.4 Le bidual E* de E est le dual de E*. 


Ex-6 On suppose que Æ est un espace vectoriel de dimension finie. Soit q 
l’application de E dans E** définie par 


(oùxeEet fe E*). 

Montrer que g est un isomorphisme de E sur E**. (Attention, ici le fait que 
E soit de dimension finie est primordial. En dimension infinie, q(Æ) n’est pas 
égal à E*). 


Ex-7 Montrer que si (f;)1<;<n est une base de E*, il existe une unique base 
(ei)icien dont (fj)1<;<n est la base duale. Aïnsi on peut parler de couples de 
bases duales. 


Ex-8 Soit Æ un espace vectoriel de dimension n. Supposons qu’il existe une 
famille (e;)1<i<n de E et une famille (f;)1<;<n de E* telles que f;(e;) = 6;,;. 
Montrer qu’alors (e;)1<;<, est une base de Æ et que donc (f;)1<,<, est sa 
base duale. 
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8.2 Exemples de formes linéaires - Applica- 
tions 


Ex-9 Montrer que pour tout 1 < à < n, l’application p,; de R? dans R qui 
à tout x = (x1,---,2,) associe x; est une forme linéaire sur E. Montrer que 
la famille (pi)1<i<n est une base de R"*. De quelle base la base (pi)1<icn est 
elle la base duale ? 


Ex-10 Soit R|X] l’espace des polynômes à coefficients réels. 

a) Montrer que si a € R, l'application 0, qui à tout polynôme P fait 
correspondre sa valeur P(a) au point a est une forme linéaire sur R[X]. 

b) Montrer que l’application 7 qui à tout polynôme P fait correspondre 
son intégrale sur le segment [—1, 1] est une forme linéaire sur R[X]. 

c) Montrer que l’application qui à tout polynôme P fait correspondre la 
valeur de sa dérivée au point a est une forme linéaire sur R[X|. 


Ex-11 (formule des trois niveaux) 

a) Notons R2[X] l’espace des polynômes de degré < 2 sur R. Rappeler la 
dimension de R2[X]. 
Notons 0, la forme linéaire qui à tout polynôme P de degré < 2 fait corres- 
pondre sa valeur P(a) au point a. Soit y un réel distinct de —1 et de 1. On 
considère les 3 formes linéaires 01,0_1,0,. Montrer qu’il existe 3 polynômes 
P:, P-1,P, qu'on calculera, vérifiant 6,(P;) = 6;,; (0 si à Æ j, 1 sinon) pour 
tout à et tout j dans {1,—1, 7}. 
En conclure que (91,6-1,0,) est une base de R,[X}* et que (P1, P 1, P,) est 
une base de R2[X]. 
Notons I la forme linéaire sur R2[X1] définie par 


Montrer qu’il existe 3 constantes À, y, v telles que 
L= UN + 1-1 + vo, 


et donc que pour tout polynôme P de degré < 2 on a 


1. P(bdt = XP(1) + wP(—1) + wP(o). 
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36 Chapitre 8 


b) On se place maintenant dans l’espace R3[X] des polynômes de degré < 3 
sur R et comme dans le a) on considère les 3 formes linéaires A, A_;, À, où 
À, est la forme linéaire qui à tout polynôme P de degré < 3 fait correspondre 
sa valeur P(a) au point a. On note J la forme sur R3[X] définie par 


Montrer que si on pose P4(X) = (X—1)(X+1)(X—7), la famille (P1, P1, P,, Pi) 
est une base de R3[X]. 

Calculer A,(P:) pour a € {1,-1,7}. 

Montrer qu’il existe une valeur et une seule de 7 telle que J(P3) = 0. On 
calculera explicitement cette valeur. 

Montrer que pour y = 70, J est dans le sous espace engendré par A1, ÀA_:,AÀ,. 
En conclure qu’il existe des constantes Ào, Lo, 2 qu’on calculera explicite- 
ment, telles que pour tout polynôme P de degré < 3 on ait 


[Pa = nPG) + Pt + Pt) 


1 


Ex-12 On se place dans l’espace R3[X] des polynômes de degré < 3 sur R. 
On considère les 4 formes linéaires #;.; où 4 et j sont dans {0, 1}, définies par 


bi3(P) = D'P(-1). 


où D est l’opérateur de dérivation. 
Montrer qu’il existe 4 polynômes F;; qu’on determinera tels que 


Pis (Prs) = dinde 


Ainsi les (P,;) et les (4;;) forment un couple de bases duales. 


Ex-13 On se place dans l’espace R.[X] des polynômes de degré < 1 sur R. 
Soient « et 5 tels que —1 < & < B < 1. On note 6, et 08 les formes linéaires 
sur R.[X] où comme dans les exercices précédents, 8,(P) = P(u). Montrer 
qu’il existe deux polynômes P,(X) et Ps(X), qu’on calculera explicitement 
tels que Ou(Pe) Tr Ou, (u, ve {a, B} ). 

En conclure que si a et 5 sont fixés, il existe des constantes À et u telles que 


TI = Na + uôg 
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où 1 désigne la forme linéaire sur R1[X] définie par 


ce qui veut dire que pour tout P € R;[X] 
I(P) = AP(a) + HP(B). 


S'inspirer des exercices précédents pour montrer qu’il existe un choiïx et un 
seul des points a et 5 pour lequel cette dernière formule est vraie pour tout 
polynôme de degré < 3. Déterminer ces points a et Ë ainsi que les constantes 
À et u qui découlent de ce choix. 


Ex-14 (polynômes de Lagrange) Sur R on considère des points 
To < Li << Tn- 


On note 6, la forme linéaire sur l’espace R,[X] des polynômes de degré < n 


définie par 8,(P) = P(u). Considérons les polynômes 

X — x)... (X — xy_1)(X — x (XX — Ty 
pa) à = ad = BE = bis) (X = 80) 
(te = xo)...(xx = zr_1)(Te = Tk41)...(Tx = de) 


Montrer que pour 4,j € {0,---,n} 
dei (L5) = dj. 


En conclure que les polynômes L£;(X) forment une base de R,[X], ayant 
(dx); pour base duale. 

Soit P(X) € R,[X]. Alors P(X) se décompose dans la base des Z4;(X) sous 
la forme 


P(X) = Na La(X). 


Calculer les coefficients a; de cette décomposition. 

En conclure qu’il existe un polynôme de degré < n et un seul qui prend aux 
n + 1 points x; des valeurs fixées y. 

Quel est ce polynôme ? (On l’appelle polynôme d’interpolation de La- 
grange vérifiant les conditions données). 


Ex-15 Soit E un espace de dimension n sur le corps fini F, à q éléments. 
Quel est le nombre d’hyperplans vectoriels distincts ? Quel est le nombre 
d’hyperplans vectoriels distincts contenant un sous espace fixé de dimension 
k<n? 
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8.3 Compléments sur la dualité 


Ex-16 Soient Æ et F deux espaces vecoriels de dimension finie et T une 
application linéaire de E dans F. On définit l’application T* du dual F* de 
F dans le dual E* de E par 


Montrer que T'+ est une application linéaire de F* dans E*. (C’est la trans- 
posée de f). 


Ex-17 Soient Æ et F deux espaces vectoriels de dimension finie et T une 
application linéaire de E dans F, q l’application de E dans E** définie par 
g(x)(f) = f(x), qi l'application de F dans F* définie de la même façon. 
Montrer que T'*(q(x)) = g1(T(x)). 


Ex-18 Soit F un espace vectoriel de dimension finie n et (f1,---, fn) une 
base de F. Soit C un sous espace de dimension k de F. Notons 6 l’application 
de F x F dans le corps de base définie par 


px, y) LL D TiVi, 
i=1 


oùx=d titi et y=d wi. Notons T l'application de F dans F* qui 
à tout y € F associe T(y) = @, de telle sorte que d,(x) = @(x, y). 
a) Montrer que T est un isomorphisme de F sur F*. En particulier pour 
tout g € F* il existe un unique y € F tel que pour tout x € F, g(x) = d(x, y). 
b) Définissons maintenant 


CT ={ye F;bj(r) =0 VreC}. 


Montrer que C € Ci. 
Montrer que si C1 € C2 alors CC Ci. 
Montrer que C+ = Ctt., 

Montrer enfin que C = C++. 


Chapitre 9 


Matrices et applications 
linéaires 


Introduction 


Nous nous plaçons dans cette partie dans des espaces vectoriels de dimension finie. Nous 
allons voir qu’une application linéaire peut être lorsque des bases sont choisies dans l’es- 
pace de départ et celui d'arrivée représentée par un tableau de coefficients qu’on appelle 
une matrice. En conséquence on peut opérer sur ces tableaux de nombres en pensant aux 
opérations correspondantes sur les applications linéaires et réciproquement. Nous expose- 


rons cette notion en nous référant en permanence à cette double interprétation. 


9.1 Applications linéaires entre espaces vec- 
toriels munis de bases 


Soient E et F deux espaces vectoriels de dimensions finies respectives m et n 
sur un même corps X. Soient (e1,e2,---,e») une base de E et (fi, f2,*--, fn) 


une base de F°. Soit enfin f une application linéaire de E dans F. Tout élément 
x de E s'écrit de manière unique sous la forme 


mm 
TX — ; Tj6; 
ÿ=1 
et donc 
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(x) = D f(e;). 


Par conséquent l’application f est connue dès que les images f(e;) sont 
connues. Réciproquement si on fixe les f(e;) la formule précédente définit 
bien une application linéaire. Se donner une application linéaire est donc 
équivalent à se donner les images des vecteurs de la base choisie. Il est com- 
mode de stocker dans un tableau rectangulaire appelé matrice les compo- 
santes des vecteurs f(e;) dans la base (fi, f2,-*-, fn). Plus précisément on 
met à la colonne j et à la ligne à la composante a;; du vecteur f(e;) sur 
le vecteur de base f;. Aïnsi la colonne j du tableau contient bien le vecteur 
(sous forme de composantes) f(e;). 


Définition 9.1.1 La matrice 


A = (a) 1<i<n 
I<j<m 
ou encore 
di, 1,2 °° im 
= G2,1 22 ‘' Am 
QAn1 An2 ‘‘' An;m 


est la matrice de l’application linéaire f dans les bases données. Cette matrice 
a n ligne et m colonnes ; on dira que sa taille est (n, m). 


Remarquons bien que si on utilise d’autres bases pour les espaces E et F'on 
change aussi la matrice de l’application linéaire f. 


Soit E un espace vectoriel de dimension finie, l’identité de E dans E,, lorsque 
on choisit la même base dans E espace de départ et dans E espace d’arrivée 
a pour matrice 


appelée matrice unité. 
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9.2 Exemples de matrices d’applications linéaires 


Ex-1 Soit E le plan vectoriel muni d’une base orthonormée. On considére 
dans £ la rotation vectorielle d’angle 0 (E est ici à la fois espace d’arrivée et 
espace de départ avec la même base). Quelle est sa matrice ? 


Ex-2 On muni R° de la base habituelle e; = (1,0,0),e> = (0,1,0),e3 = 
(0,0,1) et R° de la base f; = (1,0), f = (0,1). Soit à l’application linéaire 
de R° dans R° définie par @(x1,72,23) — (ax:,b2) (où a et b sont des 
constantes réelles). Quelle est la matrice de 6 ? 


Ex-3 Soit X,[X] l’espace des RONA de degré < n à coefficients dans 
K. On munit cet espace de la base (1,7,7+, Le -.., 2). Soit D l’opérateur 
de dérivation de K,[X] dans lui-même. Quelle est la matrice de D? 

Ex-4 Reprenons les hypothèses et les notations de l’exercice précédent. Soit 
Ty l'opérateur de translation sur K,[X] défini par T,(P)(X) = P(X +h). 
Exprimer 7, en fonction de h et de D (penser à la formule de Taylor pour 
les polynômes). 

Quelle est la matrice de T}, ? 


Ex-5 On considère le plan de la géométrie euclidienne, muni d’un repère 
orthonormé (0, e1,e2). Soit t l’application linéaire de ce plan dans lui même 
de matrice (pour la base fixée) 


in) 
(où À > 0). 


Quelle est l’image ON par t du vecteur OM de coordonées (x, y) ? 

Donner une construction géométrique de cette image. 

Soit M5 = (1,0). On construit successivement les points M, tels que OM, = = 
t(OM,-1). Montrer que les points M, sont sur une spirale logarithmique. 


Ex-6 Soit E l’espace vectoriel sur le corps {0,1} des applications affines 
de {0,1}° dans {0,1} c’est-à-dire les applications de la forme (x1,%2,x3) — 
Uo + T1 + U2X2 + Uu3æ3 dont on prend pour base les 4 applications 


(x, T2, T3) D: L 


(dis; Ts) Ti; 


(x, To, za) + To, 
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(T1, T2, T3) + T3. 


Soit F l’espace vectoriel des fonctions de {0,1}* dans {0,1} muni de la base 
(eijoci<r Où e(t1,T2,%3) = 1 si à = %1 + 2%2 + 4x3 (autrement dit si à à 
T1, T2, T3 pour développement binaire) et e;(x1,%92,x3) = 0 sinon. 

Soit Z l'application qui à f € Æ fait correspondre f € F. Quelle est la matrice 
de 7 dans les bases considérées ? 


Ex-7 Soit X3[X] l’espace des polynômes de degré < 3 à coefficients dans 


le corps K et xo, T1, to, x3 des points distincts de À. Soit T l’application de 
K3[X] dans A“ définie par 


T(P) = (P(&0), P(x1), P(&2), P(x3)). 


Quelle est la matrice de T lorsque X3[X] est muni de la base des monômes 
et K4 de sa base habituelle (1,0, 0,0), (0,1,0,0), (0,0,1,0), (0,0,0,1)? (les 
matrices de cette forme ou leurs transposées sont dites matrices de Van- 
dermonde). 

Quelle est la matrice de T lorsque K3[X] est muni de la base des polynômes 
de Lagrange (Lo, Li, Lo, La) (cf. Formes linéaires, Ex-14) ? 


Ex-8 Soit X3[X] l’espace des polynômes de degré < 3 à coefficients dans 
le corps Æ et x1,x2 des points distincts de X. Soit T l'application de X3[X] 
dans K‘ définie par 


T(P) = (P(&1), P'(x1), P(œ2), P'(æ2)). 


Quelle est la matrice de T lorsque X3[X] est muni de la base des monômes 
et K* de sa base naturelle (1,0,0,0), (0,1,0,0), (0,0,1,0), (0,0,0, 1) ? 


Ex-9 Soit u = (uo, ui, u, u3) un point fixé dans R*. On considère l’applica- 
tion T,, de R‘ dans R‘ définie par 


Ta(fo, di fa, 3) =. (go; 91; 92; 93), 


m= D, few. 


k+l=i (mod 4) 


avec 


Quelle est la matrice de T, lorsque R° est muni de sa base naturelle ? 
Remarque : Ici on peut considérer les éléments de R° comme des fonctions 
(définies par leurs valeurs) de Z/4Z dans R. Alors si f = (fo, f1, f2, fs) est 
une telle fonction, la fonction g = 7,(f) est la convolée f x u. (Les matrices 
de cette forme sont dites matrices circulantes). 


Chapitre 10 


Opérations élémentaires sur les 
matrices 


Introduction 


Il existe de nombreuses opérations intéressantes sur les tableaux matriciels. Nous présentons 
ici les opérations simples qui découlent des opérations de base sur les applications linéaires : 
l’addition, la multiplication par un scalaire, la composition, la transposition. 


10.1 Opérations sur les matrices 


Soient Æ, F,G trois espaces vectoriels sur le même corps K, munis des bases 
respectives € — (e;)iciém, Ÿ = (fijicjens 9 = (gk)i<k<p- 

Ex-1 Montrer que si a et b sont des applications linéaires de Æ dans F de 
matrices respectives (a;,;);,; et (b,,);,,, alors a+b a pour matrice (a; +bi,;)i,5. 


Définition 10.1.1 La somme de deux matrices À = (a;,;);,; et B = (bi; 
de même taille est la matrice C = (c,;);,; de même taille où c;,; = à; + bi. 


On note C = À + B. 


Ex-2 Montrer que si a est une application linéaire de Æ dans F de matrice 
(a,5)i,; et si À € K, alors Àa a pour matrice (Ai; );,;. 


Définition 10.1.2 Le produit d’une matrice À = (a;,;);; par un scalaire À 
est la matrice C = (c;);, de même taille où ci; = Aa. On note C = XA. 
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Définissons le produit d’un vecteur ligne par un vecteur colonne 

bi 

(a1, Dore An) : 

bn 
comme étant le scalaire 977, ab;. 
Ex-3 Soit b € L(E,F) une application linéaire de E dans F de matrice 
B = (b;,;) ayant donc n lignes et m colonnes, et a € £L(F, G) une application 
linéaire de F dans G de matrice À = (a;,;) ayant donc p lignes et n colonnes. 
La composée a o b est une application linéaire c de E dans G et sa matrice 
sera notée C = (c;,;). Cette matrice à p lignes et m colonnes. 


Montrer que le coefficient c;,; est obtenu en faisant le produit de la ligne 2 de 
la matrice À par la colonne j de la matrice B. 


Définition 10.1.3 Le produit de deux matrices À = (a;,;);,; de taille (p,n) 
et B = (b;,);,; de taille (n,m) est la matrice C = (c;;);,, de taille (p,m) où 
G,j — = Gi k0k,j- On note C = AB. 


Ex-4 Soit a € L(E,F) dont la matrice est À = (a;,;);; de taille (n,m). 
Considérons l’application transposée a* de F* dans E* définie par a*(b)(x) = 
p(a(x)). Montrer que la matrice (af ;);,; de a* quand F* est muni de la base 
duale de la base f et E* muni de la base duale de la base e est de taille (m,n) 
et vérifie af; = ji. 


Définition 10.1.4 La transposée de la matrice À = (a;,;) de taille (n,m), 
est la matrice A! = (a;;) de taille (m,n). 


Ex-5 Montrer que (AB)! = BfA. 


10.2 Structures des espaces de matrices 


Compte tenu de l’interprétation des matrices en termes d’applications linéaires, 
les structures de bases seront les mêmes. 

Notons Myn(K) l'ensemble des matrices ayant m lignes et n colonnes à 
coefficients dans K. Avec l’addition des matrices et la multiplication par un 
scalaire, M n(K) est muni d’une structure d’espace vectoriel sur K. 


Ex-6 Construire une base de M,,n(K). Quelle est sa dimension ? 


Opérations sur les matrices A5 


Lorsque m = n on à l’espace M,(K) des matrices carrées d'ordre n. La 
multiplication des matrice est alors une opération interne de cet ensemble et 
lui confère une structure d’anneau (non commutatif). L'unité de l’anneau est 
la matrice unité. 


Ex-7 Donner des exemples de matrices carrées de même taille À et B telles 
que AB £ BA. 


Ex-8 Montrer que si À et B sont des matrices inversibles dans M,(K) alors 
(AB)! = B71A"1. 
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Chapitre 11 


Déterminants 


Introduction 


À une famille de n vecteurs dans un espace de dimension n, ou encore à une matrice carrée 
de taille (n,n) (qui peut être considérée aussi comme la famille formée par ses n vecteurs 
colonnes), on associe un élément du corps de base appelé le determinant du système, qu’on 
peut concevoir de manière intuitive comme le ” volume algébrique” du ”parallélépipède” 
(en dimension n) construit sur les vecteurs de la famille. Ainsi ce déterminant sera 0 si et 
seulement si la famille de vecteurs est liée. 


11.1 Formes multilinéaires alternées - Déterminants 


Soit £ un espace vectoriel sur le corps K, de dimension finie n, dans lequel 
on fixe une base e = (e;)1<i<n. 


Définition 11.1.1 Une forme n-multilinéaire alternée f est une appli- 
cation de E" dans K qui est linéaire par rapport à chacune des n variables, 
et qui vérifie en outre pour tout couple 1 < i < j <n et tout (x1,---,%n) 


Fra ae ro) SES es dires 


Ainsi lorsque f est n-multilinéaire alternée, en utilisant la décomposition de 
toute permutation en produit de transpositions, on peut écrire pour toute 
permutation o de {1,---,n} 


(Go; to Et) = (—1) 0) f(x, SEA ka) 
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où e(o) est la signature de la permutation o. 


Ex-1 Montrer que l’ensemble des formes n-multilinéaires alternées est un 
sous espace vectoriel de l’ensemble des fonctions de £” dans K. Nous noterons 
An(£) cet espace vectoriel. 


Ex-2 Soit (x1,---,1,) une famille de n vecteurs de Æ, dont les composantes 
dans la base e fixée sont données par 


Soit 6 € A,(E). 


Montrer que 

(x; RAS Tn) = os (1) x 001) CCE at) o(e1, PT Én)s 

CESn 
où S, est le groupe des permutations de {1,---,n}. 
Soit À € K. Montrer que l’application 6, de E” dans K définie par 
UNETE to La) = À >. (1) x, °°" Tno(n) 
OESn 

est une forme n-multilinéaire alternée. 


En conclure que dim(A,(E)) = 1, et que toute forme multilinéaire alternée 


est définie par sa valeur prise sur (e,---,e,). 

Définition 11.1.2 On appelle déterminant dans la base e = (e1,---,e»), 
la forme n-multilinéaire alternée qui vaut 1 sur (e1,---,e»). 

Ainsi 


detftirnt) > (—1) 21,041) se Ton): 


OESn 


Le déterminant d’une matrice carrée est le determinant de la famille de ses 
vecteurs colonnes. 


Ex-3 Utiliser la formule précédente pour calculer dans le cas n = 2 


d,1 2 
det : Ë 
Q2,1 2,2 


Déterminants 49 


Ex-4 Montrer que det(A) = det(A). 
Ex-5 Montrer que det(\A) = \"det(A). 
Ex-6 Donner un exemple dans lequel det(A + B) # det(A) + det(B). 


Ex-7 Montrer que si la famille (x1,---,x,) est liée, alors det(x:,---,x,) = 0. 
En conclure que si À est une matrice carrée, et si B est obtenue à partir de 
À en ajoutant à une colonne (resp. à une ligne) une combinaison linéaire des 
autres colonnes (resp. des autres lignes) alors det(A) = det(B). 


Ex-8 Soit À = (a;;);,; une matrice carrée d’ordre n. Si 1 <i,j < n, notons 
À;,; la matrice carrée d'ordre n — 1 obtenue en supprimant de la matrice À 
la ligne 2 et la colonne 7. 

Montrer que 


det(A) = S (1) adet(Ais), 
11 
puis 


n 
det(A) = Ÿ (—1) Ya;çdet(Ai). 
j=1 
La première formule donne le développement du déterminant suivant la co- 
lonne j, la deuxième formule donne son développement suivant la ligne 1. 


Ex-9 Soit À une matrice diagonale dont les coefficients diagonaux seront 
notés À;. Quel est le déterminant de À ? En particulier quel est le déterminant 
de la matrice unité ? 


Ex-10 Montrer que det(AB) = det(A)det(B). 
En conclure que si À est inversible det(A) À 0 et que det(A7!) = 1/det(A). 


Ex-11 Montrer que det(A) # 0 est équivalent à À inversible. 


Ex-12 Soit Æ un espace vectoriel de dimension n muni d’une base e. Soit 
(fj)1<j<m une famille finie de vecteurs de £. On appelle À la matrice de 
taille (n,m) dont les vecteurs colonnes sont les vecteurs f; décomposés dans 
la base e. Montrer que le rang du système de vecteurs (f)1<;<m est r si et 
seulement s’il existe une matrice carrée B extraite de À de taille (r,r) de 
déterminant non nul, et toute matrice carrée de taille (r + 1,7 + 1) extraite 
de À est de déterminant nul. 
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11.2 Calcul de déterminants 


Les problèmes de calcul de déterminants se classent dans deux catégories 
distinctes. D’une part le calcul de type formel pour les déterminants qui ont 
une forme particulière, d'autre part le calcul numérique de déterminants sans 
forme algébrique particulière. Le cas du calcul formel de déterminants n’est 
pas un cas d’école. De nombreux problèmes concrets conduisent en effet à 
s’intéresser à des déterminants ayant une forme particulière, et pour lesquels 
on arrive à obtenir une formule close. 

Le calcul numérique de déterminants ” quelconques” se fait dès que l’ordre 
dépasse trois ou quatre par des méthodes numériques adaptées (par exemple 
pivot de Gauss). 


Ex-13 Calculer le déterminant 
1 


1 QT: 1 


Ex-15 Calculer le déterminant 


x y T+Y 
y T+Y Z 
T+Y TZ y 


Ex-16 calculer le déterminant d’ordre n tridiagonal en fonction des pa- 
ramètres à,b,c 
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a b 0 0 ... 0 0 

ce a b 0 ... 0 0 

0 € a b ... 0 0 
1e 

b 

C à 


On pourra établir une relation de récurrence liant D,, à D,,_1 et D, — 2. 


Ex-17 (Déterminants de Vandermonde) On considère le déterminant 


RC EL 
V(2, 21;°°'; Zu) = 


On remarquera que ce déterminant est un polynôme en les n + 1 variables 
20, **", Zn. Quel est le degré de ce polynôme ? 

Montrer que ce polynôme à pour tout à < j le monôme (2; — z;) en facteur. 
Calculer ce déterminant. 


Ex-18 Calculer le déterminant 


1 1 1 .. 1 
il (821 2 tt An-1 

2 2 2 
Di |t 4 RO 
n—1 n—1 n—1 
1 à Q An_1 

où 
2ikr 
AK — de 


Ex-19 Calculer le déterminant circulant d’ordre n 


Go Qi *‘"" An-1 
QAn—1 ao °t"  An-2 
(24 & ... (#7) 


(on pourra multiplier à droite ce déterminant par le déterminant de l’exercice 
précédent). 
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Ex-20 Calculer le déterminant d’ordre 2n +1 


1 cos(xo) sin(xo) --- cos(nto)  sin(nxo) 
1 cos(xi) sin(xi) --- cos(nx:)  sin(nx) 
1 Cos(ton) Sin(ton) **+ COS(NTon) SiN(Ntan) 


Ex-21 On veut calculer le déterminant D, d'ordre n dont le coefficient 
d'indice (4, j) est ne (on supposera que a; + b; # 0). 
Montrer dans un premier temps que D, est de la forme 


[ai — a;)(bi — b;) 


D, = Cy mn 
[1 (ai + b;) 


où C, est une constante. 
Montrer que 


ai +b1 a1+bn 
lim a, D, = ; 
An —00 dr. she 1 
An—1+01 An—1+0n 
puis que 
1 1 0 
a1+b1 a1+0n—1 
lim lim a, Dy = 
bn —00 An —}00 he, us 1 
An 1+01 An—1+bn-1 
En conclure que 
| . _ AnD 
lim lim ==], 


bn — 00 An —00 


puis que pour tout n, Cn=1. 


Chapitre 12 


Applications linéaires en 
dimension finie 


Introduction 


Nous reprenons ici l’étude des applications linéaires d’un espace de dimension finie dans 
un autre espace de dimension finie. Nous regarderons alors ce qu’apporte la notion de base 
à cette étude. 


12.1 Bases et applications linéaires 


Soit E et F deux espaces vectoriels sur le même corps K, de dimensions res- 
pectives m et n. Soit e = (e1,---,em) et f = (f1,---, fn) des bases respectives 
de E et de F. 6 est une application linéaire de E dans F. 

Ex-1 Montrer que Im() = [@(e1),---,p(em)|. 


Définition 12.1.1 On appelle rang d’une application linéaire, la dimension 
de son image. 


Ex-2 Soit G un sous espace vectoriel de E. Montrer qu’il existe une appli- 
cation linéaire 0 de E dans F telle que Ker(0) = G. 


Ex-3 Soit G un sous espace vectoriel de E. Dans Æ on définit la relation 
xRy si et seulement six -yeG. 
Montrer que R est une relation d'équivalence. 
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Montrer que sur l’ensemble quotient on définit une structure d’espace vecto- 
riel en posant 


cl(x) + cl(y) = cl(x + y), 
Àcl(x) = d(Ax), 


où cl(x) désigne la classe d'équivalence de x. 


Définition 12.1.2 L'espace vectoriel définit précédemment est appelé es- 
pace quotient de E par G, on le note E/G. 


Ex-4 Soit D un supplémentaire de G dans E. 
E=GeD. 


Si x et y sont deux éléments de E alors x et y s’écrivent de manière unique 
sous la forme x = gx + dy et y = y + dy AVEC 9x, Jy € G dx; dy € D. 
Montrer que si x — y € G alors dy = dy. 

On peut alors définir l’application p de E/G dans D par 


p(cl(x)) = ds. 


Montrer que p est une application linéaire. 

Montrer que p est bijective. 

Quelle est la dimension de E/G? 

On remarque qu’il est en un certain sens équivalent de travailler avec un 
supplémentaire de G ou avec E/G. 


Ex-5 (factorisation des applications linéaires) Montrer qu’il existe une 
application injective ® unique de E/Ker(®) dans F, telle que 9 = & o s où 
s est la surjection naturelle de E sur E/Ker(). 

51 de plus @ est surjective alors ® est bijective. 


Ex-6 Soit D un supplémentaire de Ker(@) dans Æ. Montrer que la restriction 
®: de 6 à D est surjective. Comparer ®&; à l’application ®& de l’exercice 
précédent. 


Définition 12.1.3 Soit e! = (e!,---,e!,) une autre base de E. La matrice P 
de l’application identité de E muni de la base e! sur E muni de la base e est 
la matrice de passage de la base e à la base e’. Ainsi les colonnes de P 
sont les coordonnées des vecteurs de la nouvelle base e' dans l’ancienne base 
e. 
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Ex-7 Montrer que P est inversible. 


Ex-8 Soit e = (e!,---,e,,) une autre base de E et f! = (f,,---,f!) une 
autre base de F. Notons P la matrice de passage de la base e à la base e’ et 
Q la matrice de passage de la base f à la base f’. Si À est la matrice de 6 
pour les bases e et f, montrer que la matrice À’ de @ pour les bases e/ et f’ 
vérifie 

A=Q "AP 
Ainsi, dans le cas où F = E et où f = e, f! = e' on obtient pour l’endomor- 
phisme 

A! = P-TAP. 


Ex-9 Etudier le cas où F = E*, f = e* (base duale de E) et f! = e*. 
Si P est la matrice de passage de e à e’, quelle est la matrice de passage de 
e* à e* ? 
Montrer que 
A' = P'AP. 


Ex-10 Soit E un espace vectoriel de dimension finie n, e = (e1,---,em) une 
base fixée dans Æ, f une application linéaire de E dans E (endomorphisme ou 
opérateur) de matrice À dans les bases données. Montrer que À est inversible 
(f inversible) si et seulement si det(A) Z 0. 


Ex-11 Expression de l’inverse d’une matrice. 

Soit À = (a;,;);,; une matrice carrée d'ordre n. Si 1 < ,j < n, notons À;,; la 
matrice carrée d'ordre n — 1 obtenue en supprimant de la matrice À la ligne 
à et la colonne 7. 

On appelle cofacteur de a;,; le nombre (—1)*/det(A; ;). 

Montrer que 


où C;,, est le cofacteur de a;;. 


12.2 Exemples 


Ex-12 Soit R,[X] l’espace des polynômes de degré < n. Soit m < n et 
T0, * *  Xm des réels distincts. Soit F le sous espace des multiples du polynôme 
N(X) = (X-x0)---(X-x,,). Quelle est la dimension du quotient R,[X|/F ? 
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Déterminer un isomorphisme de ce quotient sur R,,[X|. 


Ex-13 Dans le plan de la géométrie euclidienne on considère une base or- 
thonormée (e:,e2). Soit À la rotation vectorielle d’angle 0. On considère la 
base (f1, f2) où f1 = Role) et f2 = Ro(e2). Quelle est la matrice de passage 
de la base (e:,e2) à la base (fi, f2) ? On considère l’application linéaire de 


matrice | 
or —pisin(6:) ) 
pisin(:)  picos(®:) 


dans la base (e:,e2). Quelle est la matrice de cette application linéaire dans 


la base (fi, f2) ? 
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Systèmes linéaires 


Introduction 


La résolution d’un système linéaire de m équations à n inconnues (pouvant prendre leurs 
valeurs dans un corps donné) s’interprête en terme d’algèbre linéaire comme la recherche 
d’un vecteur dans un espace de dimension n dont l’image par une application linéaire 
donnée est un vecteur donné dans un espace de dimension m. Nous verrons à ce propos 
plusieurs problèmes distincts : existence d’une solution, unicité, détermination de la ou 
des solutions, algorithmes effectifs de calcul. 


13.1 Problèmes linéaires en dimension finie 


Un système linéaire de m équations à n inconnues est un système du type 


Gaiti ++ @inln = bi 
Q2,1T1 a un dnln — bo 
Gm,1T1 Ré Zmnln  — bn 


où les inconnues x; sont à chercher dans un corps K, les coefficients a;,; du 
système sont dans À, et les données b; des seconds membres sont aussi dans 
K. 


Pour un tel système les problèmes posés sont les suivants : 
a) Le système possède-t-il au moins une solution ? 
b) S’il existe une solution est elle unique ? 
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c) S'il existe des solutions peut on les écrire suivant une formule en 
fonction des coefficients et des seconds membres ? 

d) Comment calculer numériquement une solution lorsqu’ on se donne 
des valeurs numériques explicites pour les coefficients et les seconds membres ? 
En particulier ce calcul peut il être réalisé pratiquement à partir de formules 
closes ? 


La première chose que nous allons faire est d’interpréter de diverses manières 
en termes d’algèbre linéaire le problème posé. Nous verrons que ces diverses 
façons de voir sont toutes intéressantes par les divers éclairages qu’elles 
apportent au problème. 


La première présentation est celle sous forme de système que nous avons 
donnée. Partant de là, considérons l’application linéaire L de K” dans K” 
qui à tout x = (x1,%2,---,%n) fait correspondre y = (y1, Yo, -- +, Ym) OÙ Y% = 
= &,j2,. Notons b = (b:1,b2,---,b,,) le vecteur de K7* constitué par les 
seconds membres du système. Dans ces conditions on cherche les solutions de 
L(x) = b, c’est-à-dire les vecteurs x de X” dont l’image par L est un vecteur 
donné b de K”. Sous cette forme il est clair que la résolution du système 
linéaire passe par l’étude de l’application linéaire L. D'ailleurs, si on munit 
K7 et K7 de leurs bases naturelles, on peut écrire le système sous forme 
matricielle 


Az = b 
ou encore 
1 192 ‘°" Ain T1 b: 
2,1 22 ‘° An To bo 
Gm,1 Am2 ‘°° Amn Ta bn 


Ainsi que nous l’avons déjà remarqué dans le chapitre sur les formes linéaires, 
on peut considérer que l’étude d’une application linéaire de K7 dans X7? peut 
se faire par l’étude de m formes linéaires sur 7. Plus précisément, notons 
L; (1 <i< m) la forme linéaire sur K" définie par 


L;(x) — ÿ di,jli: 
j=1 


Dans cette optique la recherche des solutions du système se ramène à la 
recherche des solutions communes aux m équations 
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Comme les solutions de chaque équation définissent un hyperplan affine 
(cas L; Æ 0), ou l’espace tout entier (cas L; = 0,b; = 0), ou encore l’en- 
semble vide (cas L; = 0,b; 0), on voit bien la signification géométrique 
de l’ensemble des solutions : c’est une intersection d’objets du type précédent. 


Remarquons encore que si on appelle À;, 45,---, À, les vecteurs colonnes 
(donc éléments de X7*) de la matrice À, le problème posé se ramène aussi 
à trouver les coefficients x; d’une combinaison linéaire de ces vecteurs qui 


vaille b : - 
DETTE 
i=1 


Enfin compte tenu de l’interprétation du problème en terme d’algèbre linéaire 
on peut penser à l’exprimer de manière qui semble plus générale : 


Soient Æ et F deux espaces vectoriels de dimensions finies respectives n et 
m sur un corps K, et L une application linéaire de Æ dans F'. Etant donné 
be F existe-t-il x € E tel que L(x) = b? 


En fait en fixant une base dans Æ et une base dans F, on se ramène à un 
système du type défini précédemment, si bien que cette façon de poser le 
problème n’est pas réellement plus générale. Cependant elle est intéressante 
car elle permet de s’affranchir de l’aspect analytique qu’imposent des bases 
fixées pour une expression plus intrinsèque. A titre d'exemple citons le problème 
d’interpolation suivant : considérons Uo,---,umn_1 des points distincts de 
R. Soit L l’application linéaire de R,;_:[X] (espace des polynômes sur R 
de degré < n—1) dans R°” qui à tout P € R,_:[X] fait correspondre 
(P(uo),---, P(um_1)). Etant donné b = (bo,---,bm-1) un élément de R”, 
trouver P tel que L(P) = b (ou encore tel que pour tout 0 < à < m — 


Remarquons encore que le problème sous sa forme intrinsèque peut se présenter 
à l’aide de formes linéaires sous la forme suivante : soit £ un espace vecto- 
riel de dimension n sur K; soient ZL:,---,L,, des formes linéaires sur Æ et 
b1,---,b,, des éléments de Æ ; trouver les vecteurs x de Æ tels que pour tout 
1<i<m on ait Li(x) = b;. À titre d'exemple le problème d’interpolation 
cité au dessus s’exprime facilement sous cette forme. 

En conclusion disons que nous utiliserons pour chaque problème linéaire à 


résoudre l’une ou l’autre de ces formes, en choisissant bien entendu la plus 
commode compte tenu du contexte. 
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Ex-1 On reprend les notations de l'introduction et on suppose m = n. 
Montrer que les conditions suivantes sont équivalentes : 


1) le problème linéaire L(x) = b à une solution et une seule pour tout 
be KT”; 

2) Lest bijective: 

3) Ker(L) = {0}: 

4) det(A) # 0; 

5) il existe b € K°7° tel que le problème L(x) = b ait une solution et une 
seule. 
Donner des exemples et des contre-exemples. 


Définition 13.1.1 Lorsque, dans le cas où m = n, les conditions équivalentes 
précédentes sont réalisées, on dit que le système linéaire considéré est un 
système de Cramer. Remarquons que le fait d’être un système de Cramer 
ne dépend pas du second membre du système. 


Ex-2 Montrer que si m £ n il n’est pas possible que le système ait une 
solution et une seule ( c’est-à-dire que soit il n’a pas de solution, soit il à 
plusieurs solutions). Donner des exemples de chaque cas. 


Dans le cas général l’analyse du problème n’est pas très difficile à faire en 
regardant le noyau et l’image de l’application linéaire L. 


Ex-3 Si le système a une solution x = (x1,---,x,) l’ensemble de toutes les 
solutions est x + Ker(L). 

Le système a une solution si b € Im(L). En particulier le système homogène 
associé (L(x) = 0) a toujours une solution ( 0 est toujours solution). 

Si le rang de L est m (dimension de l’espace d'arrivée F) le système a 
toujours une solution, sinon il existe des cas où il n’en a pas. 

Si le noyau de L est réduit au vecteur nul il y a au plus une solution. 


Les divers cas précédents peuvent éventuellement être détectés par des valeurs 
de déterminants. 


Ex-4 Nous cherchons à résoudre L(x) = b où L est une application linéaire de 
l’espace vectoriel E de dimension n dans l’espace vectoriel F de dimension m 
(espaces dans lesquels on suppose que des bases ont été fixées). Nous notons 


A = (ai) 1<jen là matrice de L dans les bases fixées. Soit r le rang de L 
1<j<m 
(dimension de l’image). On sait alors qu’il existe un sous déterminant d’ordre 
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r non nul extrait de la matrice À (on supposera pour fixer les idées que c’est 
celui obtenu pour 1 <i<r,1<3j<r)et que tout sous déterminant d’ordre 
r + 1 extrait de À est nul. Montrer que le système admet une solution si et 
seulement si pour tout 1<s<m—7r 


1,1 AU dir b: 
2,1 Lt d2,r bo 
0. 
Qr,1 RÉ Arr b, 
Qr+s,1 CU, Qr+sr bis 


Ex-5 Dans le cas d’un système de cramer dont la matrice (carrée) est 


A = (@,j) jen dont les seconds membres sont b.,---,b, l'unique solution 
1£ÿ£n 
(T1, ---, Th," , Zn) S’écrit 
Qi, ‘‘" @1,k—1 bi Q1,k+1 ‘°' in 
_— ni ‘‘  An,k-1 bn An,k+1 ‘°° Ann 
RS 
di ‘‘°" in 
ni ‘‘ Ann 


Les formules précédentes, dites formule de Cramer ne sont pas employées 
pour la résolution explicite numérique des systèmes (sauf éventuellement en 
petite dimension). Parmi les méthodes explicites de résolution, la méthode 
du pivot de Gauss (ou encore méthode par élimination) est très simple. 
Divers algorithmes plus ou moins élaborés réalisent cette élimination. Nous 
proposons ici l'algorithme suivant dit algorithme de Crout. 


Ex-6 Soit À une matrice carrée d’ordre n dont le déterminant est non nul 
À = (@,j)i<i,jen 
On cherche deux matrices carrées d’ordre n 
L = (ohiéiiés avec ;,j — 0sj>1 


U=(Bishisisen avec bi; = 0 si j <i 
telles que À = LU. 
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Combien y a-t-il de coefficients à déterminer ? 
Combien y a-t-il d'équations qui permettent de les déterminer ? 
Ceci explique qu’on puisse en outre imposer 


Bii=1 1<i<n. 


Puisqu'on sait que les coefficients B;; valent 1 , on ne les stocke pas. On 
va donc pouvoir mettre les deux matrices L et U dans le même tableau T 
carré d'ordre n (L matrice triangulaire inférieure et U matrice triangulaire 
supérieure). 

Pour chaque couple (i, j) vérifiant 1 < 4, 7j < n on dispose donc de la relation 

inf(i,5) 
QikÜkG = Gi 
k=1 

étape 1 


1) Calculer les a;1 en fonction des à; 1. 

2) Permutation des lignes afin que parmi tous les a; 1 ainsi calculés, celui de 
plus grande valeur absolue se trouve en position (1,1). Par abus de notation 
nous appellerons toujours ces nombres à;1. Montrer qu'après permutation 
dQi,1 F 0. 

3) Calculer les B,; pour j > 1 en fonction des a:,; et de @1. 

On peut itérer ce calcul : supposons calculées les p — 1 premières lignes et les 
p — 1 premières colonnes. 


étape p 


1) Calcul de la p° colonne (à partir du rang p) c’est à dire des éléments 
a, pour à > p. Montrer que à;, se calcule en fonction de a;, et d'éléments 
déjà calculés. 

2) Permutation des lignes pour mettre le a;, (où à > p) de plus grande 
valeur absolue en position (p,p). On montrera alors que le @,, ainsi trouvé 
est non nul. 

3) Calcul de la p° ligne à partir du rang p + 1, c’est à dire des éléments 
Bp,j pour j > p. Montrer que B,,; s'exprime en fonction de a,,; et d'éléments 
déjà calculés. 


Remarque 


Lorsqu'on utilise cette transformation de la matrice À sous la forme LU 
pour résoudre le système AX = B on est amené à tenir compte des diverses 
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permutations faites sur les lignes en faisant des permutations analogues sur 
les composantes de B. on a alors ensuite à résoudre LUX = B, soit encore 
LY = B puis UX = Y. 


Ex-7 Montrer que le système L(x) = b admet une solution si et seulement 
si pour toute solution u de L*(u) = 0 on a u(b) = 0. 


13.2 Résolution explicite de systèmes linéaires 


Ex-8 Résoudre 


27%: — To T3 —= 4 
3T1 +4%9 —2%3 =: ‘11 
3T —2%9 +43 = Ii 


Ex-9 Quand le système suivant possède-t-il une solution unique ? 


AX1 +bt2 +Ct3 +dxza = 0 
bx, —at2 +dxz3 —cxa —= 0 
Ct) —dto —ax3 +bxa —= 0 
dti <+Cta —bxz3 —-axa —= 0 
Ex-10 Résoudre le système 
Ti —2To +X3 + = 1 
Æi —2T2o +X3 —X4 — —1l 
T1 —2To +Z3 +5%4 = 5 


Ex-11 Résoudre le système 


Ti +To —3%3 = —] 
27%: + To —2%3 —= 1 
Ti +Lo +T3 —= 3 
Ty +2To —3%3 — ll 


Ex-12 Soit à résoudre un système de Cramer du type 


d1T1 +- + Ginln b: 


QAn,1T1 DEEE Annln br 
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Notons pou 1<k<n 


la solution du système 
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d,1T1 Hier dinln => "| 
Qg-1121 + +@k-1nEn — Ù 
GkA1T1 + °°" + Ann = 1 
dpt Fe oeneta =: 0 
Anti + °°: + Anntn = 0 
Montrer que la solution cherchée est 
DSi +--- + bSn. 
On pourra montrer que (S:,---,S,) est une base de X7. Quelle est sa base 


duale ? 
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Valeurs propres, vecteurs 
propres 


Introduction 


Nous étudions ici les notions très importantes de vecteurs propres et valeurs propres. 
Nous traiterons d’abord le cas simple de la diagonalisation des applications linéaires. Nous 
reprendrons ensuite l’étude de la décomposition spectrale de manière plus complète, en se 
plaçant toutefois dans le cas où le corps de base est algébriquement clos. 


14.1 Une première étude des valeurs propres 
et vecteurs propres 


Définition 14.1.1 Soit E un espace vectoriel de dimension finie n sur C. 
Soit T un opérateur de E. Un vecteur propre de T est un vecteur non 
nul e tel que T(e) = Àe où À € C. Un tel nombre À est appelé une valeur 
propre de T. 


Ex-1 Montrer que pour que À soit valeur propre de T il faut et il suffit que 
det(T — AT) = 0. En conclure qu’il existe au moins une valeur propre et au 
plus n valeurs propres. 


Définition 14.1.2 Le polynôme en À det(T — AI) est appelé polynôme 
caractéristique de T.. 
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Ex-2 Soit V, l’ensemble des vecteurs propres associés à la valeur propre À. 
Montrer que S\ = VA U {0} est un sous espace vectoriel de Æ£. Cet espace est 
appelé sous espace propre associé à À. 


Ex-3 Montrer que si 7 vecteurs de Æ sont des vecteurs propres associés 
à r valeurs propres distinctes de T, alors ces r vecteurs sont linéairement 
indépendants. 


Ex-4 Montrer que si le polynôme caractéristique de T n’admet pas de 
racines multiples alors il existe une base de Æ formée de vecteurs propres 
de T (attention le corps de base est € ou un corps algébriquement clos). 


Ex-5 Supposons qu'il existe une base de Æ formée de vecteurs propres de T 
(par exemple comme nous l’avons vu si toutes les valeurs propres de T' sont 
distinctes). Quelle est la matrice de T dans cette base ? 


On à ainsi dans l’exercice précédent diagonalisé la matrice de T (c’est- 
à-dire trouvé une base de Æ dans laquelle la matrice de l’opérateur T est 
diagonale). Plusieurs problèmes se posent alors. Que se passe-t-il dans le cas 
où le polynôme caractéristique a des racines multiples ? Peut on toujours 
réduire la matrice de Ten choisissant une bonne base ? Nous allons pour 
répondre à ces questions développer le paragraphe suivant. Il restera à voir 
ce qu'il se passe quand le corps n’est pas algébriquement clos, et que le 
polynôme caractéristique n’a pas autant de racines qu’il le faudrait. Nous 
ne développerons pas cette question. Disons simplement que dans ce cas, on 
peut grossir” le corps de base en se plaçant dans sa clôture algébrique. Ceci 
demande de définir la ’complexification” d’un espace vectoriel, ce que nous 
ne ferons pas ici. 


14.2 Décomposition spectrale des opérateurs 


Soit Ÿ un espace vectoriel de dimension finie n sur le corps © des nombres 
complexes. Soit T un endomorphisme de # et A(T) la sous algèbre de 
l'algèbre L(X) des endomorphismes de # engendrée par l’opérateur T. 


14.2.1 L’algèbre engendrée par un opérateur. 


Il est facile de voir que l’algèbre A(T) est constituée par les opérateurs P(T) 
où les éléments P sont les polynômes de l’espace C[X]. On remarquera que 
A(T) est une algèbre commutative. 
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Soit ® l’application de C[X] dans A(T) qui à un polynôme P fait corres- 
pondre l'opérateur P(T). Cette application ® est visiblement surjective. Elle 
ne peut pas être injective car A(T') est de dimension finie alors que C[X] 
est de dimension infinie. Le noyau Ker(®) de ® n’est donc pas réduit à {0}. 
Ce noyau est un idéal de C[X]. On en déduit que Ker(®) est engendré par 
un polynôme M normalisé de degré minimal unique (en effet C[X] est un 
anneau principal). Remarquons que ce polynôme n’est pas nul (le noyau n’est 
pas {0}) et n’est pas 1. Donc le degré de M est supérieur ou égal à 1. On 
peut donc énoncer le théorème suivant : 


Théorème 14.2.1 1! existe un polynôme normalisé et un seul M tel que 
tout polynôme P vérifiant P(T) = 0 soit un multiple de M. Ce polynôme de 
degré minimal parmi ceux qui sont annulés par T est de degré supérieur ou 
égal à 1. 


Définition 14.2.1 Le polynôme M du théorème précédent est appelé po- 
lynôme minimal de l’opérateur T.. 


Puisque le corps de base est algébriquement clos il est clair que le polynôme 
minimal se décompose en produit de facteurs de la forme (X — À). Ceci nous 
incite à étudier les opérateurs de la forme (T — AT). 


14.2.2 Les opérateurs (T — X)f. 


Pour tout 4 > 0 et tout À € C définissons 


NŸ = Ker(T — \1)* 


Il est clair que la suite de sous espaces (WŸ);>0 est croissante (pour l’inclu- 
sion). Comme % est de dimension finie, cette suite ne peut être strictement 
croissante, il existe donc un plus petit entier (À) appelé index de À tel que 


NN ce A 


Théorème 14.2.2 Si (À) est l'index du nombre complexe À alors pour tout 
entier p > 0 


Ne _ NRQUE 
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Démonstration Montrons par récurrence que pour tout p > 0 N AP = 


N, APTE, On sait que le résultat est vrai pour p=0. Supposons le vrai pour 


un entier p > 0 et montrons alors qu’il est vrai pour p +1. Six € NOR 
alors (T — M)(x) € NX *P*T et par hypothèse de récurrence (T — AT)(x) € 


NEC Donc x € | AAA et par suite NEC = NP 


Notons n1 la dimension du sous espace W ; a), Compte tenu de la croissance 
stricte de la suite (WF), pour 0 < k < v(À), on a l'inégalité v (À) < n1 < n. 


14.2.3 Spectre, index et polynôme minimal. 


Définition 14.2.2 Le spectre de l’opérateur T, noté o(T), est l’ensemble des 
nombres complexes À tels que T — ÀT ne soit pas inversible. 


Remarque : 

— Nous sommes ici en dimension finie, donc il est équivalent dans ce cas de 
dire que le spectre est constitué des À tels que T' — ÀZ ne soit pas injectif, 
c’est à dire tels que Wi Æ {0}. 

— Si À n’est pas dans le spectre de T alors l’index (À) est nul. Par contre 
s’il est dans le spectre, son index est supérieur ou égal à 1. 


Définition 14.2.3 Lorsque À est une valeur propre (ie. v(À) > 1), le sous 
espace N° est appelé sous espace propre de T et ses éléments non nuls sont 
les vecteurs propres de T.. 

Le sous espace NS est appelé sous espace spectral associé à À. 


Théorème 14.2.3 Le spectre de T est un ensemble fini et le polynôme mi- 
nimal M de T s'exprime sous la forme 


M(X)= ][ (X-2® 


ÀA€o(T) 


Démonstration Soit À un polynôme non nul tel que R(T') = 0. Supposons 
que À ait une racine À qui ne soit pas dans le spectre. Alors R(X) = (X — 
ÀA)R(X). Ainsi R(T) = R;(T)(T — ÀÏ). Et comme T — ÀJ est inversible 
et R(T) = 0 on en conclut que R1(T) = 0. Si R à une racine À qui est 
dans le spectre, et qui a une multiplicité (À) + p avec p > 1 alors R(T) = 
(T — A)OPR(T). Donc R(T) = 0 implique que pour tout r R?(T)(x) 
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appartient au noyau de (T — A1)/@)+? c’est à dire au noyau de (T — AI)", 
On a donc aussi en posant R3(T) = (T — M)NR(T) l'égalité Ra(T) = 0. 
La conclusion de cette première partie est qu’il existe un polynôme annulé 
par T, donc multiple de M, qui a pour racines des éléments À du spectre avec 
une multiplicité < (À). Soit donc maintenant P un tel polynôme annulé par 
T'et soit À un élément du spectre de T. Alors il existe un vecteur y non 
nul tel que T(y) = Ày. Si bien que P(T)(y) = P(A)y, et comme P(T)=0, 
on obtient P(ÀA) = 0. (On vient de montrer en passant que le spectre est un 
ensemble fini). Supposons maintenant que la multiplicité a de cette racine 
À soit strictement inférieure à (A). Alors il existe x1 % 0 tel que (7 — 
AD) (x1) = 0 et y1 = (T — AI)%(x1) # 0. P(X) peut s’écrire sous la forme 
P(X) = (X — A)Q(X) avec Q(À) £ 0. Comme T(y1) = Àyi on peut écrire 
P(T)(x1) = Q(T)(y) = Q(Ü)y £ 0 ce qui contredit le fait que P(T) = 0. 


14.2.4 Décomposition de l’espace en sous espaces spec- 
traux. 


Théorème 14.2.4 Chaque sous espace NA est stable par l'opérateur T.. 


Démonstration T(T — A)/® = (T — AI) AT donc si x € Ker(T — A)" 
alors T(x) € Ker(T — A)". 


Lemme 14.2.1 Si À £ y alors NY AN: = {0}. 


Démonstration En effet soit x un élément supposé non nul de l’intersection 
et à < v(À) le plus grand entier tel que z = (7 — A1)*(x) # 0. Donc 
T(z) = Àz. Si bien que (T — p1)"W(2) = (À — p}"Wz Z 0. D'autre part 
(T2 a (2) = (7 = OT = MD(a) soit (T WI} W(e) = (T - 
AI)(T — p1)"W (x) = 0. D'où la contradiction. 


Les projecteurs associés au polynôme minimal. 


Introduisons les polynômes Q;\(X) = ns Les polynômes (Q\)1e0(7) 
sont premiers dans leur ensemble, d’après le théorème de Bezout il existe 


donc des polynômes À, tels que 


ÿ_ AQ;=1 


À€o(T) 
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Posons pour tout À du spectre 
S(X) = AUX)Q(X). 


Les polynômes S, vérifient les propriétés suivantes : 


Théorème 14.2.5 


1. La somme des S\ vaut 1. 


2. Le produit S\S, lorsque À Æ u, est un multiple du polynôme minimal 
M. 


3. Le produit de S\ par (X — Xe) est un multiple du polynôme minimal 
M. 


4. S\(À) = 1 et pour tout autre élément 1 Æ À du spectre S\(u) = 0. 


Démonstration Les propriétés 1) 2) et 3) proviennent immédiatement des 
définitions des polynômes S,. La propriété 4) provient de 1) et de la définition 
de 5. 


Notons alors pour tout À du spectre J\ = S\(T). Les opérateurs 1, vérifient : 


Théorème 14.2.6 


1 Dxeoçr) Lx = 1: 

2. 3 À À u alors T1, = 0. 
STE 

4 (T — ADO, = 0. 
ETL=LT 


Démonstration 1) 2) et 4) sont des conséquences immédiates du théorème 
précédent sur les polynômes S,. la propriété 3) est une conséquence claire de 
1) et 2). La propriété 5) provient du fait que /) est un polynôme en T. 


Remarque : Les opérateurs 7, sont des projecteurs. 
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Décomposition de l’espace. 

Notons #, = J,(#). 

Théorème 14.2.7 L'opérateur T laisse stable X,. 
Démonstration Ceci provient de TI, = IT. 


Théorème 14.2.8 On a pour tous les éléments À et 1 distincts, dans le 
spectre les relations suivantes : 


14) = {0} 
et pour tout x € X; 
HG=T 


Démonstration Ceci est une conséquence immédiate du fait que les opérateurs 
I, sont des projecteurs qui vérifient Z,1, = 0. 


Théorème 14.2.9 L'espace se décompose sous la forme 


4= x 
À€o(T) 
Démonstration Du fait que 3 '\égçr) 1 = 1'il est clair que 4 = 5 é(r) A. 
Si D_xco(T) xx = 0 avec x1 € #,. Pour tout y dans le spectre on peut écrire 
LD xeotr) TA) = Zn D'où x, = 0 et ceci pour tout y. 
Théorème 14.2.10 
(A 


Démonstration Puisque (1 —A1)"%)1, = 0 on conclut que #\ C NA. Soit 
alors x € NY CE D'après le résultat précédent x s’écrit de manière unique sous 
la forme x = et z1 Où æ, est dans #, (et par suite aussi dans NI), 
Donc 


0=(T-M)O(z)= NO (T- M) (x). 
HEo(T) 
Mais les (T — A7)"%)(x,) sont eux même dans #, à cause de la stabilité par 
T de ces sous espaces. Donc en vertu de l’unicité de la décomposition de 0 
on en conclut que pour tout y dans le spectre (T — A1)"W(x,) = 0. Mais si 
À ju comme x, est dans NE il ne peut être aussi dans N° que s’il est 
nul (lemme 5.1). Par suite x = x} et donc x € 4. 
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Définition 14.2.4 Les projecteurs 1, sont appelés les projecteurs spectraux 
de T.. 


Remarque : 

— Désormais pour étudier l’opérateur T on pourra l’étudier sur chacun des 
sous espaces spectraux NY N, Nous verrons par la suite des exemples de 
cette stratégie. 

— La somme des dimensions n\ des sous espaces NY Q) est n. La somme des 


index est donc < n. 


Définition 14.2.5 Un opérateur est dit diagonal si Æ est la somme directe 
de ses sous espaces propres. 


Remarque : Ceci se produit si et seulement si les index des valeurs propres 
sont tous égaux à 1. Ce qui veut encore dire que dimKer(T — AI) = n, pour 
toute valeur propre À. Il est aisé de voir qu’un opérateur Test diagonal si et 
seulement s’il peut se mettre sous forme diagonale dans une base bien choisie. 


Théorème 14.2.11 Tout opérateur T s'écrit de manière unique sous la 
forme T = D + N où D est un opérateur diagonal et où N est un opérateur 
nilpotent qui commute avec D. 


Démonstration Définissons l’opérateur D par D(x) = À er) Ata où & = 
D_xcocr) ta. Soit alors N = T—D. Il est clair que la décomposition T = D+N 
est bien du type indiqué. Pour l’unicité supposons qu’on ait pu décomposer 
T' sous la forme indiquée. D étant diagonal on peut trouver des sous espaces 
E;,E2,---,E,% tels que 4 = Œ. E; et que la restriction D, de D à E; soit 
Xl. Six € E, alors T(x) = À;x + N(x) ou encore N(x) = (T — À;l)(x); 
N commutant avec D il laisse stable les sous espaces propres E; de D, donc 
désormais tout se passe dans E;, en particulier N'(x) = (T —-À;1)"(x). Comme 
N est nilpotent il existe z, tel que (T — À;1)"#(x) = 0. En prenant x # 0 on 
conclut que Ker(T — À;1)"# Æ {0}, donc que À; € o(T). Par suite il existe 
À € o(T) tel que E; C N\ a (et deux E; distincts sont dans des NY Le 


distincts). Par suite les E; sont les NY ® et on retombe sur la décomposition 
dejà trouvée. 
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14.2.5 Polynôme caractéristique 


Définition 14.2.6 On appelle polynôme caractéristique de T le polynôme 
A(X) = det(T — XI). 


Compte tenu de l'interprétation du determinant, les racines de À sont les 
valeurs pour lesquelles T — XJ n’est pas injectif, donc sont les valeurs propres 
de T. Il suffit de décomposer T sur les sous espaces NY de. pour voir que 


A(X)= L[.4=2" 


À€o(T) 
En particulier comme (À) < n, À est un multiple de M donc 


Théorème 14.2.12 {Cayley Hamilton) Le polynôme caractéristique est an- 
nulé par T.. 


14.2.6 Appendice : Une autre approche des projec- 
teurs spectraux 


Fonctions analytiques d’ opérateurs 


La construction que nous avons donnée des projecteurs spectraux utilise le 
théorème de Bezout sur les polynômes. On va se placer ici dans une classe 
de fonctions plus large afin d’obtenir une partition de l’unité plus directe. La 
classe que nous allons utiliser est celle des fonctions analytiques. Le premier 
pas est donc de définir l’image d’un opérateur par une fonction analytique et 
d’ étudier quelques propriétés des opérateurs ainsi obtenus. La définition en 
vue utilise l’interpolation de Lagrange Sylvester. 


Notons #(T°) la classe des fonctions d’une variable complexe à valeurs com- 
plexes qui sont analytiques dans un ouvert contenant le spectre de T (cet 
ouvert pouvant varier suivant la fonction considérée dans la classe). 


Si f € H(T) on sait qu’ il existe un polynôme P tel que pour tout À € a(T'), 
et tout m < v(À) —1 
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Il est clair que la définition est indépendante du polynôme d’ interpolation 
P (cf. théorème 4.1). 
Les résultats suivants découlent immédiatement des propriétés des polynômes. 


Théorème 14.2.13 Si f et g sont dans H(T) et si a et B sont deux nombres 
complexes, alors 


1. af + Bg E HIT) et (af + Bg)(T) = af(T) + Bg(T). 
2. fa € HT) et fa(T) = f(D)a(D). 

CR DIN Ne An 0e US DD 

4 


. f(T) = 0 si et seulement si f)(X) = 0 pour tout À € o(T) et tout 
0O<m<v(À)-1 


Il découle de 2) que les opfateurs f(T) et g(T°) commutent. 


Les projecteurs spectraux 


Définissons pour tout À € a(T°) un disque ouvert D) de centre À de telle sorte 
que ces disques soient deux à deux disjoints. Pour tout À € o(T') on définit 
alors sur l’ouvert réunion des disques ainsi construits la fonction f\ qui vaut 
1 sur D, et O0 ailleurs. Il est clair que ces fonctions sont analytiques. On note 
alors Æ; l'opérateur f\(T). 
Théorème 14.2.14 Pour tout À € o(T) on a 

E, — TI, 


c’est-à-dire que les E\ sont les projecteurs spectrau. 


14.3 Exemples 


Ex-1 Trouver le polynôme caractéristique, les valeurs propres, les vecteurs 
propres et diagonaliser (si possible) les matrices suivantes 


5 2 
A= | 101 
Le 1200 
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sé 4. 0 0 
LT OÙ 0 

20 .-1l. «6; 12 
LG. "6 “M 


2 1 
(42) 
Ex-2 Trouver le polynôme caractéristique, les valeurs propres, les sous es- 


paces propres des matrices suivantes. Sont elles diagonalisables ? Quels sont 
les polynômes minimaux des endomorphismes correspondants ? 


1 1 1 
L = = 
ll 4, = 


1 
1 
Fe 1 
ES ES 


Ex-3 On appelle bloc de Jordan toute matrice de la forme 
À 1! 
À 1 0 


NE 
0 À 


ou éventuellement une matrice (À) de taille 1. 


Montrer que tout opérateur linéaire d’un espace vectoriel complexe Æ de 
dimension finie admet une matrice constituée de blocs diagonaux de Jordan 
et de 0 ailleurs. 

Quelle forme simple peut on avoir si Æ est un espace vectoriel sur R ? 


Ex-4 Mettre sous forme de Jordan la matrice B de l'exercice 2. 

Ex-5 Soit T un opérateur linéaire d’un espace vectoriel complexe Æ de 
dimension finie. Soit À une valeur propre de T et supposons que le sous 
espace propre (ker(T — ÀT)) associé à À soit de dimension 1. Montrer que 
l'index v() est égal à la dimension n, du sous espace spectral WŸ Lee 
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Chapitre 15 


Interventions de l’algèbre 
linéaire 


Introduction 


Nous donnons dans ce chapitre, quelques exemples d'utilisation de l’algèbre linéaire dans 
divers domaines. Nous avons choisi comme illustration la théorie des signaux discrets finis 
et une partie de l’analyse numérique élémentaire. 


15.1 Fonctions sur un groupe abélien fini, théorie 
du signal fini 


On fixe un entier n > 2. Notons G le groupe Z/nZ et F l’espace des fonctions 
définies sur G à valeurs dans le corps des nombres complexes C. 
15.1.1 L’espace F 
À tout élément a de G on associe la fonction e, définie par 
= 1 si zx =a, 
FC era 
1) Montrer que la famille (e4)o<a<n-1 est une base de F. Aïnsi F est un 
espace vectoriel sur €, de dimension n. 


2) Soit f € F. Quelles sont les composantes de f sur la base (ec )o<a<n-1. 
3) Décrire la base duale de la base (e,)o<a<n-1: 


IT 
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15.1.2 Les caractères de G 


A tout élément u de G on associe la fonction x, définie par 


2imuv 


YU) = Er + 


1) Montrer que les fonctions x, vérifient 


Xe +9) = Xu(r)xu (y), (15.1) 
Xu(0) = 1. (15.2) 


2) Montrer qu’on obtient ainsi toutes les fonctions complexes définies sur 
G vérifiant (1) et (2). 


15.1.3 Produit scalaire hermitien sur F 


1) Montrer qu’on définit sur F un produit scalaire hermitien en posant 


< f,g >= F(u)g(u). 


ueG 
2) Montrer que la base (es )o<a<n-1 est orthonormée. 


3) Calculer pour tout u et tout v dand G le produit scalaire < X4, Xy >. 
En conclure que la famille (X4)o<u<n-1 est une base de F, et que cette base 
est orthogonale. 


4) Calculer en fonction des valeurs prises par f, les composantes de f dans 
la base (X4)o<u<n-1: 


15.1.4 Transformation de Fourier 
Pour toute fonction f € F, on définit la fonction Fe F, par 


Fo) = fer. 


ueG 


1) Montrer que 


ft) = 2 fo 


veG 


En conclure que 
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») 


es 


(u) = nf(-u). 

2) Comparer < f,g > et < f,9 >. Comparer la norme de f avec celle de 
Te 

3) Notons F l’opérateur de F qui à f associe f. Quelle est la norme de 
l’opérateur F. 


4) Calculer lorsque Jf =e, et lorsque f = x4. 


15.1.5 Matrices associées aux objets précédents 


1) Quelle est la matrice de F si F est muni de la base (es )o<a<n-1 ? 


2) Quelle est la matrice de F si F est muni de la base (Xu)o<u<n-1 ? 
3) Calculer Det(F). 


4) Calculer la matrice de passage P ainsi que son inverse, de la base 
(eos à la base (Xa)o£uéne ts 


15.1.6 Fonctions de F et polynômes formels 


À tout fonction f € F on associe le polynôme P; € C[X] de degré < n —1 
défini par 


Montrer que 


et que 


En conclure un algorithme pour trouver le polynôme P d’interpolation de 

Lagrange, de degré < n — 1, qui prend des valeurs données aux n points 
2iru 

en. 

15.1.7 Convolution et filtres stationnaires 


On définit la convolée f: x f de deux fonctions de F par 
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fi x f(x = NA +y)p(- y). 


yEeG 


1) Montrer que 


fix f(x) = fe x fix = Ÿ fi(u) 


U+V—=X 


2) Calculer fx hp et fi. 


3) Etudier P,;,,, en fonction de P;, et P;,. 

4) Pour tout t € G soit T; l'opérateur de translation qui à la fonction f 
fait correspondre la fonction 7;(f) telle que 7,(f)(x) = f(x +t). Soit H un 
opérateur qui commute avec 7, pour tout t. Montrer que les x, forment une 
base de vecteurs propres de H. Quelles sont les valeurs propres associées ?. 
Montrer qu’il existe une fonction h € F telle que pour tout f € F, on ait 


H(f)=hxf. 


15.2 Analyse numérique : Interpolation 


15.2.1 Introduction 


L’interpolation est un sujet très vaste lié aux questions d’approximation 
des fonctions. Très grossièrement il s’agit de trouver dans une classe fixée de 
fonctions (par exemple les fonctions polynômiales) un élément réalisant un 
certain nombre de contraintes. Souvent ces contraintes sont liées à la donnée 
d’une fonction f qu’on cherche à approcher par un procédé d’interpolation 
(par exemple la fonction cherchée doit prendre la même valeur que f en 
des points donnés). On se trouve alors confronté à plusieurs problèmes de 
natures différentes. Tout d’abord un problème algébrique, celui de trouver 
le ou les éléments de la classe choisie qui réalise les contraintes. Ensuite un 
problème d’approximation qui consiste lorsqu’on est parti d’une fonction f à 
mesurer la qualité de l’approximation théorique obtenue. Enfin un problème 
algorithmique, celui de déterminer un algorithme performant qui permette de 
calculer facilement et de manière aussi exacte que possible la ou les solutions. 


Dans un premier temps nous partirons d’un exemple important : l’interpola- 
tion de Lagrange. 
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15.2.2 Interpolation de Lagrange 


Soient To, T1, …, Zn des nombres complexes distincts et Yo, y1, …, y, des nombres 
complexes. Il s’agit de trouver un polynôme P(X) vérifiant P(x;) = y pour 
toutes les valeurs de k comprises entre 0 et n. 


L’aspect algébrique 


Existence de solutions Notons C[X] l’espace des polynômes à coefficients 
complexes et C,;[X] le sous espace des polynômes à coefficients complexes de 
degré inférieur ou égal à n. Considérons alors l'application linéaire T de C[X] 
dans C"*! qui à un polynôme P(X) fait correspondre (P(xo), P(x1), .…, P(xn)). 
On voit que le noyau Ker(T') de l’application Test l’espace constitué des 
multiples du polynôme N(X) = (X — xo)(X — x1)...(X — x,) et qu’on peut 
écrire 


CIX] = CA[XT D Ker(T). 


Ceci nous montre que la restriction de T à C,[X] est une bijection de C,[X] 
SHC 


Le résultat obtenu est donc le suivant 


Théorème 15.2.1 Pour tout élément (yo, Y1, .…, Yn) de C"*? il existe un po- 
lynôme P(X) de degré < n et un seul (polynôme d'’interpolation de Lagrange) 
tel que P(xx) = y (0 <k <n). Tout polynôme Q(X) (de degré quelconque) 
vérifiant aussi Q{(xx) = yx (0 < k < n) s'écrit sous la forme 


Q(X) = P(X) + K(AÏN(X): 


Ecriture sous la forme naturelle (base des monômes) Si on cherche à 
trouver explicitement le polynôme de Lagrange écrit sous la forme habituelle 
P(X) = ao + a X +. +a,X" on est amené à résoudre le système de n +1 
équations en les n + 1 inconnues &o, …, An 


Go + MiTo +. + AnEg — Yo 
QT Eee 'alT = Hi 


0 Fin + ee Ont = Un 
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Ce système est un système de Vandermonde dont on sait bien sûr qu’il ad- 
met une solution unique. Nous fournirons ultérieurement un algorithme pour 
résoudre ce système plus rapidement que ne le ferait une méthode classique 
comme la méthode du pivot par exemple. 


Ecriture sous forme de Lagrange Ainsi le calcul des coefficients a; du 
polynôme cherché se ramène à la résolution d’un système de Vandermonde. 
Sans être très compliqué ceci n’est cependant pas immédiat. Or il peut se faire 
qu’on n’ait pas absolument besoin des coefficients a; et qu’on puisse se sa- 
tisfaire d'exprimer le polynôme d’interpolation dans une base mieux adaptée 
que la classique base des monômes. Dans cet ordre d’idée, introduisons les 
n + 1 polynômes ZL;(X) (où 0 < k < n) qui vérifient 


Le(xy) = 0 GA) 
D’après l’étude qui précède le polynôme L,(X) existe et est unique. On vérifie 
aisément que L4(X) s’écrit explicitement sous la forme 


(X —= ul CX —= zk-1)(X — Tri (X — ds) 
(tk == 20): (Ex = De-1)(Ee = Lila (TE == Lu) ‘ 


Ces polynômes forment une base de C,;[X] et le polynôme d’interpolation 
s'exprime dans cette base 


Lx(X) = 


P(X) = due La(X). 
k=—0 
Il est intéressant de remarquer que le polynôme L;(X) s'écrit aussi 


N(X) 
(X — zx) N'(xx) 
Net) ee) 


Le(X) = 


Ecriture sous la forme de Newton L’inconvénient des polynômes L;(X) 
est que chacun d’eux fait intervenir tous les points d’interpolation et donc si 
on rajoute un nouveau point tout calcul fait à partir des polynômes L;(X) 
doit être entièrement refait. 

Prenons alors les polynômes NX) = (X — %0)(X — z1)..(X — xyx_1) où 
0<k<n (avec N = 1). Ces polynômes forment aussi une base de C, [X] et 
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le polynôme d’interpolation se décompose sur cette base sous la forme de 
Newton 


Le problème est alors de calculer les coefficients b;. Pour ce faire définissons 
les différences divisées successives des valeurs y; par rapport aux points x; 


[vo] = Yo 


[y1, res x] E" [vo, Gus Yk_1] 


[vo, V1,..., Uk] — 
Tr — X0 


Théorème 15.2.2 Les coefficients de la décomposition du polynôme d’ in- 
terpolation de Lagrange de degré n dans la base de Newton sont donnés par 


by = [Vo, V1, …, Yxl 
où 
CLEAN: 


Preuve. La formule à démontrer est clairement vraie si on a n = 0. Sup- 
posons la formule vraie pour les polynômes de degré n — 1 interpolant en n 
points. Soit alors P,_1(X) le polynôme d’interpolation de Lagrange associé 
aux points Zo, T1, .…., T1 et aux valeurs Yo, V1, …., Yn-1 Qn-1(X) le polynôme 
d’interpolation de Lagrange associé aux points 21,%9,...,7, et aux valeurs 
V1, V2, +, Un Et 


P(X) = D bNe(X) 


le polynôme d’interpolation de Lagrange associé aux points Z0,%1,..….,2, et 
aux valeurs Yo, Y1, …, Yn. Il est facile de voir que 


Pa-1(X) = 5 be Ne(X) 


si bien qu’il reste simplement en vertu de l’hypothèse de récurrence à établir 
la formule pour le coefficient b,. 
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Pour cela définissons 


Ê (E =) AA) (= Tn)Pn1(X) 


FX) = 
(X) PRES 
On vérifie que : 
PE) — Yi 
pour 0 << n. Donc : 
P(X)=P(A): 


En égalant les coeflicients du terme de degré n dans l’expression des po- 
lynômes P, et P, on obtient la relation cherchée. 


En comparant l’expression de Newton du polynôme P, d’interpolation en 
k +1 points avec celle de Lagrange on trouve en regardant les coefficients des 
termes de degré k 


k 
V5 
[Vo; 1; .… Yxl = AT LR 
2, Ni) 


j=0  k+i 
Il est clair dans cette représentation que le rajout d’un nouveau point ne 
fait que rajouter un nouveau terme au polynôme, les autres termes restant 
identiques. 


L’aspect approximation 


Le théorème de division des fonctions différentiables Soit f une 
fonction réelle définie sur R de classe C?*!, où p est un entier naturel. On 
suppose que f s’annule en un point a de R. Posons 


g(x) = f'(a + (x — a)u)du 


alors g(x) est l’unique fonction continue telle que 


f(x) = (x — a)g(x). 
De plus d’après le théorème de dérivation sous le signe intégrale on voit que 
la fonction g(x) est de classe C? et que pour tout 0 <q <p 


1 
gx) = / TF0 (a+ (x — a)u)du 
0 
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et par suite 


Fa 


1 
gta) (x)| =“ DE 1 "Ptela,r] 


Majoration de l’erreur Soit f une fonction de classe C"*!, P le polynôme 
d’interpolation de Lagrange qui prend les mêmes valeurs que f aux point 
To, Z1,-.., En et À un intervalle compact contenant x, xo, %1,…,2,. Appliquons 
alors le théorème précédent à f(x) — P(x). On obtient 


f(x) — P(x) = (x — t0)go(x) 


avec 


COPA. 
HSE 


(ne pas oublier que P+ (x) = 0) 


super | F0 (0) 


puis 


go(x) = (x — t1)g1(x) 


avec 


n— 1 n 
gi (a)| < =sumer|g" () 


et ainsi de suite. Si bien que 


LG = PO le = ne 1) — an )lsuprer1070) 


(n +1 


L’aspect algorithmique 


Le calcul des différences divisées L’algorithme des différences divisées 
est très simple. Il utilise l'écriture du polynôme d’interpolation sous la forme 
de Newton. Les coefficients sont alors calculés par la formule de récurrence 
établie précédemment 


[vo] = Yo 


a, …, xl — [Yo, …, Vk-1] 


[vo. Vis... Ve] + 
Tr — Lo 
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si bien que le calcul se fait conformément à la figure 15.1. 
Le nombre d’opérations à effectuer est en O(n?). 


[vol 
D 
[Vo, vi] 
a = 
[y] [Vo 1: y] 
[va we] . [Vo: ÿ1: Ya, al 
[ul af [ua Vo, 
[u, ya] 4 
[vs] 


F1G. 15.1 — Le calcul des différences divisées 


Résolution d’un système de Vandermonde Nous avons vu que le cal- 
cul effectif des coefficients du polynôme d’interpolation de Lagrange dans la 
base naturelle des monômes passe par la résolution d’un système de Van- 
dermonde. Voici un algorithme qui permet de résoudre un tel système. Cet 
algorithme est basé en fait sur l’algorithme de Hôrner pour l’évaluation de 
polynômes. Il est plus rapide que les algorithmes directs de résolution des 
systèmes linéaires généraux comme la méthode du pivot de Gauss ou la 
méthode de Householder qui sont en O(n*) alors que nous obtenons ici un 
algorithme en O(n°). 
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Soit N un entier > 2. Etant donné x = (x1,%2,...,xn) un N-uplet de réels 
deux à deux distincts on note B la matrice 


1 1 Ne 1 
A Lo TN 
2 2 2 
B — TT To TN 
_1 N-1 1 
Ti To TN 


On cherche à résoudre le système 


BW =Q 
où Q est la matrice colonne constituée des seconds membres q1,q2,---,qn du 
système et où W est la matrice colonne constituée des inconnues w:1, Wo, :--, Wn 


du système. 


Pour tout entier j vérifiant 1 < j < N on pose 


à ZT —X 

P;(x) =[T=——> 
a: n 
nAi 


P; est donc un polynôme en x de degré N — 1 qui peut s’écrire 


N 
Pc) = D Art 
k=1 


En effectuant le produit de la matrice A = (4;4);x par la matrice B on 
constate que 
AB = (P;(tx))5x 


ce qui prouve que À est l’inverse de B. 
On peut alors écrire que W=AQ. On obtient ainsi les formules 


N 
k=1 


Nous allons dans la suite mettre en place une méthode de résolution qui 
calcule les coefficients des polynômes P;, donc qui calcule l'inverse de la 
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matrice B.Pour calculer les coefficients de P; on sera amené à calculer les 


coefficients de 
N 


N;() = [[(& - 2) 
nai 
et aussi le dénominateur intervenant dans la formule qui définit P;, c’est à 
dire le nombre N;(x;). 
Posons 
P(x) = (x — x1)(x — x2)...(x — zn) 


P(x) est donc un polynôme de degré N qui s’écrit sous la forme : 
P(x) = 2" +ena +... + ox + 


Montrons tout d’abord comment si on connaît les coefficients c; on peut 
calculer les coefficients du polynôme 


Pour cela posons 
N,;(x) —= ba iii box + b. 


On vérifie immédiatement sur l’expression de N;(x) que le coefficient du 
terme de plus haut degré est 1. En remarquant que P(x) = N;(x)(x—x;) on 
établit la formule 

by_1 = Cr + x 0x. 


Si bien que 
bn = 1 
Di = Cr + Tijd 
Connaissant les coefficients de N,; il est alors facile de calculer le dénominateur 


N;(x;) intervenant dans la définition de P,. 
En effet posons ty = bn = 1 et définissons pour tout k < N 


ti = Lite + bx-1 


On constate alors que 1 = N;(x;), le calcul proposé pour N;(x;) n'étant rien 
d’ autre que l’algorithme de Horner. 
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Il reste maintenant à calculer les coefficients c; de P. 
Pour tout entier k vérifiant 1 < 4 < N on définit 


Qu) = (x — r1)(x — 2)... (x — xx) 
et on écrit Q, sous la forme 
Quiz) = 2° + OT NT Hills + QRRATT +... + 


Il est facile de voir sur l'expression de Q:(x) = x — 1 que 11 = —%1. 
De la formule 


Qk(t) = Qx1(x)(x — Ti). 


découlent pour k = 2,3,..., N les formules 
Akk — Ak—1,k-1 — Tk 


Op,j = Ak-1,5-1 — TkOk-1,5 Ï =k—-1,..,2 
ce qui achève l’algorithme. 


On peut voir que le nombre d'opérations à faire dans cet algorithme est en 
O(N?), la partie la plus coûteuse étant le calcul des coefficients cz. 


Un cas particulier : les points d’interpolation sont les racines n° de 
l’unité 


Interpolation de Lagrange et transformée de Fourier discrète On 
trouvera dans [5] des précisions sur la transformation de Fourier discrète. 
Rappelons que si a = (ao, &1, .…, än-1) est une suite finie de n nombres com- 
plexes on définit la transformée de Fourier & de la suite a comme étant la 
suite @ = (Go, 1, .…., dn-1) OÙ 


La transformation inverse se calcule facilement (cf. [5]) par 


n—1 

Fe 2imuv 

Ay = ) Qye  . 
v=0 
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On remarque que 
1 


(GARE, PE 
À un coefficient près le même algorithme permettra de calculer la transformée 
de Fourier et son inverse. 


Notons 
PA(X) = (ü SE Re x") 
On voit alors que | 
au = Per ). 

Cette dernière remarque met en évidence un aspect très important de la 
transformée de Fourier discrète : l’aspect interpolation . En effet il est facile 
de trouver grâce à ce que nous avons vu le polynôme d’interpolation de 
Lagrange qui prend les valeurs a, aux points zy = et : C’est le polynôme 
P;(X) dont les coefficients sont donnés par la transformée de Fourier discrète 
de la suite a = (ao, &1, …, aän-1) des valeurs prises aux points d’interpolation 


Le calcul explicite : transformée de Fourier rapide Il existe divers 
façons proches les une des autres de calculer une transformée de Fourier 
discrète. Toutes ces variantes sont des algorithmes de transformée de Fourier 
rapides (FFT). Nous nous placerons ici dans le cas où le nombre d’éléments 
de la suite à transformer est n = 2”. Pour tout r > 0 et tout 0 < k < 27 —1 
posons 


2ikr 
T 


WÈ =e-7 


Remarquons que 
Wr = (Win) = (WT 


par exemple 
(Ws) = (Wa) = Wi 


On rappelle que si 
a = (ap, @1,.…., Gom_1) 
et si ; 
PA) = 5e (ao +aX ++ am X?"-1) 
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alors 
GPO) 


Pour tout polynôme 


P(X) — Po + pi nn PhrmXe 


notons 
T—1 
Po(X) = po + paX +. + par2X7 
et 
r—1 
P(X) = pi + paX +. + par iX7 
alors 


P(X) = P(X?) + X P.(X?) 
ce qui donne si0 <k<2-1—1] 
P(We) = Po(Wa-1) + War P(Wa-1) 


et 
P(WE2 7) = P(W) — Wé P(Wg-). 


Ces dernières formules vont nous donner un algorithme pour calculer les 
valeurs de la transformée de Fourier. 


Remarquons tout d’abord que si on a tabulé les valeurs de W*, alors on 
dispose aussi des valeurs de W£ pour tout r < m. 


Pratique du calcul. Le coefficient = n’interviendra qu’à la fin. Pour cela 
au lieu de calculer avec le polynôme P, nous calculerons avec P = nP, = 
Go +. + GE ER NE 

L'exemple m = 3 est suffisamment instructif pour décrire l’algorithme. Re- 
marquons que 


P500(X) _ ao; Pooi(X) — @a, Po1o(X) = a, Poi(X) — 6 
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Q5 a3 a7 


O O 
L 
O O 
L 
O O 
| 
O4 O4 O4 O4 Ok O\ Ok OX 
ao Qi a2 a3 A4 Q5 A6 a7 


F1G. 15.2 - La FFT sur 8 points 
P:00(X) =. &, Pro (X) = a5, Po(X) = a3, Pr(X) — Q7. 


On commence donc à faire une permutation o des éléments 
Go, 1, A2, @3, 4, @5, A6, Q7 
pour les mettre dans l’ordre 
A0, 4, 2, 6, G5, 43, A7. 


Ceci se fait facilement en remarquant qu’à chaque indice supposé écrit en 
binaire on fait correspondre l'indice obtenu en écrivant les bits dans l’ordre 
inverse. Ainsi l’indice 4 = 100 est transformé en 1 = 001. la suite du calcul 
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de la transformée de Fourier se fait en trois étapes indiquées par la figure 
15.2 et à la fin on divise les coefficient obtenus par 8. 


Appelons M, M2, M$ les matrices 


L Es, cs 20 0 0 0 D 

de 2 60, 210 LA: KG 0, 0 

de A0: Gi, Gb O7, KO A0 VU 

1 D Ont CO MO ee Ù 
HE dE AO 9-40 EL 1 A E 
Dr 0: 0". 0: 1° «ef. 0" © 

de Aÿ:; 40% 01 10 D LE 

D: L'Oie 0, 406.20 07 SE, = 

de Ù- E 06 0,0: 4 0 

do Eee 0 << 07 400 4 

M 10h 20: ‘@ .04 201 "D 

2 D «cf: Os 0 0 0" 
Me D Hobc 0e ue É. 201, VE - 
D AÙs 0 20 29 dE -0 À 

Dr 07 ds 0: CE 0e = "D 

De 0 OÙ O7 LE 0 I 

Loco 20 © E C0: “nv 0 

Qe 2Le #0 30 00%: St ' 0 

D. 0 «° 20 507 «0 EE "4 

3 Date nee D CO Gi 
Me Le 0 07 sets À A (D 
D «fe 0 “0 «0: =T 0 Ù 

D 0 Le 20 207 cc ET D 

De 00e CON C0 01 0: ET 


Si, S2, S3 les matrices diagonales définies par 
Si = Diag(1,W?,1,W2,1,W9,1,WS) 
S2 = Diag(1,1,W?,W},1,1,W?,W) 
S3 = Diag(1,1,1,1,W£,W2,W2, Wè) 
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et enfin © la matrice de la permutation ” reverse bit” o. 
Dans ces conditions la matrice F3 de la transformation de Fourier sur 8 points 
s'écrit 

hs 5 1 


Ceci se généralise facilement pour n = 2”. Le nombre d’opérations à effectuer 
pour calculer cette transformation est de l’ordre de nlog(n). 


15.2.3 Le problème général de l’interpolation 


Interprétons de manière un peu plus algébrique le problème d’interpolation 
de Lagrange. En particulier notons g; la forme linéaire sur C, [X] qui à tout 
polynôme Q(X) fait correspondre Q(x;). On cherche alors un élément P(X) 
de C,[X] qui réalise g;(P) = y; pour tout i. Remarquons que les g; forment 
une base du dual de C,[X] et que cette base est la base duale de la base 
constituée par les L,;. 


Nous poserons le problème général de l’interpolation en ces termes : 


Soit E un espace vectoriel de dimension n, E* le dual de E, (g;); une base 
de E* et (y); des nombres. Trouver un élément P de E tel que g;(P) = 
pour tout 1 << n. 


Il est clair que si (e;); est la base de Æ dont (g;); est la base duale, alors 


i=1 


Nous allons voir par la suite quelques exemples qui entrent dans ce cadre 
général. 


15.2.4 Quelques exemples importants 


Interpolation d’Hermite 


Soient zo < æ1 et Yo, Y1, Yo, Y\ des nombres réels. On cherche un polynôme P 
de degré < 3 tel que 


P(xo) — % 
P(t) = y 
P'(&o) = % 
P'{&) = y 
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Notons Æ l’espace vectoriel de dimension 4 des polynômes de degré < 3. 
Définissons les 4 formes linéaires sur Æ 


00,x0 (©) P(xo) 
dom (P) — P(x) 
O1x(P) — Pt) 
Om(P) = P'{x) 


La question posée est donc de résoudre le problème linéaire suivant : trouver 
P tel que 


0x0 (P) — Vo 
Obs (P) — Yi 
Ô1,&0(P) Ce Yo 


Pour cela cherchons si on peut trouver 4 polynômes H6 36; How; 1,50? H1,m 
tels que 

1sii=ketj=l 

0 sinon L 


die (Hyx) — { 


Si on arrive à trouver ces polynômes cela prouvera à la fois que ces polynômes 
forment une base de Æ, que les formes linéaires introduites forment une base 
de E*, qui est la base duale de la base polynômiale trouvée. 


Un simple calcul nous permet de trouver effectivement ces polynômes (ce 
sont en fait des solutions dans des cas particuliers bien choisis du problème 


posé) : 


CE 
Hot) = (= a} (Bn 20 29) 
Rae Sa 
FE ne 


Il ressort de toutes ces considérations que le problème à une solution unique 
(polynôme d’interpolation de Hermite) donnée par 


96 Chapitre 15 


P(x) _ Yo Ho,x0 (T) SUR V1 Ho, (x) + Vo H1,20 (€) sd Yi Hi (x). 


Si nous sommes partis d’une fonction f de classe C{ sur un intervalle compact 

T contenant les points xo,x1, et si nous appelons P le polynôme d’interpola- 
tion de Hermite associé aux points x, x, et aux valeurs f(xo), f(x1), f'(xo), f'(x1), 
alors comme dans l’exemple de l’interpolation de Lagrange, l’application 
répétée du théorème de division des fonctions différentiables nous donne 
l’approximation 


| 


f(x) - PGI < 2 - 20) (x — a)" sup FMOIE 


& 


Interpolation par les splines cubiques 


Soit À un partage d’un segment [a, b|, 
aA—= TL <a <-- <Xn =D 
et Y = (y1,%,:--,Yn) une suite de n nombres réels. 


Définition 15.2.1 Une fonction f définie sur [a, b] est une fonction spline 
cubique relativement au partage À et à Y si les conditions suivantes sont 
satisfaites : 
1) f est de classe C?{a, b] ; 
2) f coïncide avec un polynôme de degré 3 sur chaque intervalle [x;,x;11] ; 
HI) = ui pour = 12: 
Nous noterons Say l’ensemble de ces fonctions splines. 


Nous nous proposons de déterminer des fonctions splines cubiques. 

Soit donc f € S4y et posons M; = f"(x;). Sur l'intervalle [x,,x;:1] la 
dérivée seconde de f est de degré 1 et on a sur cet intervalle 

T—T; 


Lit 4 
x) = M — + M, —<— 
Tr. — 7: J — 

Tj+1 — Lj Tj+1 — TL; 


En intégrant deux fois, 


fa)= Ge 2, Mn ba) : 00) 


6 (thx) 6 (xyn—t;) 
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où Q(x) est un polynôme de degré 1 que l’on peut présenter sous la forme 
Q(x) = a(x;n — x) + B(x — x;). 
Evaluons à et $ en utilisant f(x,;) = y; et f(x;11) = y;+1. On trouve 


= ( MO) 


6 Tjt1 — Tj) 
(ru = 35) 1 
fs M SE QE — 
B (us: j+1 6 œyu x)’ 
si bien que si on pose 
hui = Tin —%; 


on obtient pour x E[x;,x;41], 


Mbnres) a Mab=s) 


Re — — 
fG) 6 hj+i 6 hu 
hihi (in — &) kih\ (x — x) 
y; — M, ) J Es (u — M... — ) | 
( 3 J 6 Rjn j+1 j+1 6 Hs 


Nous allons exploiter maintenant les conditions sur la dérivée première. Pour 
2<j<n—-lona 


h: h: Fe). 
Re in Re y EE RE 
xs) = M. 2 LM. + 2179 
J (x;) j-1 6 + J 3 sa h; 
si bien que 
LE PE œ h;+ lin y, F hj41 Mie it Vi Yi Yi 
6 hjy h; 


On dispose donc de n — 2 équations linéaires pour déterminer les n inconnues 
M, M, ---,M,. Il convient donc si on espère avoir une solution unique de 
donner deux conditions supplémentaires. Pour la suite du calcul nous impo- 
serons donc les valeurs de la dérivée de f en x, et en x, c’est-à-dire 


Fm) — y! 
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On a alors à 
V2 U1 ! 
2M, + M = — = 
1 2 ma ( 1e à) 
et à 
Un — Yn-1 
M, 2M= | 1 
1 Pa (a, An ) 
Si on pose 
À = I 
En 
h; : 
Ne PU TRE 2<j<n—-l 
Hi = 1—); 1<j<n 
et 


— 6 | y2=y: ! 
He (se -0) 


ù = (u-e) 
alors le système à résoudre est 

2 À 0 O0 -.. -.. 0 M; bi 
Mn De 22 0 er b M bo 
0 u3 2 }3 - : | 
De 2 Fes : : = 
LDC An_2 0 

: : lt 2 ÀAn_1 î : 
Loto mate 0 a 2 M, bé, 


Ce système est tridiagonal, à diagonale strictement dominante; il possède 
une solution et une seule. On peut employer pour un tel système une sim- 
plification de la méthode du pivot, qui dans ce cas particulier s'exécute en 
temps linéaire. 


Les fonctions splines cubiques minimisent approximativement l’énergie de 
flexion d’une tige. En effet si la fonction g(x) représente l’équation d’une tige 
parfaitement élastique son énergie de flexion est 
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' (9"(0) ‘dt 
CO 


et donc si la dérivée de g(x) est petite on peut approcher cette énergie de 


flexion par 
b 2 
Î (e'œ) dt. 


a 
Nous allons voir que les fonctions splines cubiques minimisent cette dernière 
expression. Sur C?[a, b] on définit un semi-produit scalaire par 


b 
re / f"(Og" dt 


qui donne la semi-norme 


= (furet) | 


Soit WA le sous espace de C?[a, b] constitué des fonctions nulles aux points 
x; et telles que f’(a) = f'(b) = 0. 


Lemme 15.2.1 Les espaces NA et SA,y sont orthogonaux et d’intersection 
réduite à {0}. 


Lemme 15.2.2 Toute fonction f € C?[a, b] s'écrit d’une façon unique sous 
la forme f=s+hoùse Say et he Mi. 


Théorème 15.2.3 Parmi toutes les fonctions f € C?[a, b] qui vérifient f(x;) = 
y; et les conditions f'(a) = y! et f'(b) = y, la fonction qui minimise 
l’énergie de flexion est la spline cubique (l'élément de Say) qui vérifie en 
outre J'(a) = yet. f'(b) =:v;: 


On pourrait aussi montrer le résultat suivant 


Théorème 15.2.4 Parmi toutes les fonctions f € C?[a, b] qui vérifient f(x;) = 
y; la fonction qui minimise l’énergie de flexion est la spline cubique (l'élément 
de Sa,y) qui vérifie en outre f"(a) = f"(b) = 0. 
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Interpolation par des polynômes trigonométriques 


Etant donnés 2n + 1 nombres réels distincts 
—T LT <Lo << Toni LT 


et 2n + 1 nombres complexes w1,---,%Y2:1 on cherche un polynôme trigo- 


nométrique 
nm 


PÉ)= D (aycos(pz) + bsin(pz)) 


p=0 
tel que P(x:) —= Vi. 


On pourrait étudier directement ce problème, mais ici nous allons plutôt le 
relier au problème de l’interpolation de Lagrange que nous avons déjà étudié. 
Pour cela posons pour tout -n <k <n 


1 7. 
Ck — (ox — iSigne(k)bw) 


alors 
n 
Pi = > cet 
k=-n 
et aussi 
2n 
HE 24 = | Ck ere 
k=0 


En posant X = e° on obtient 
2n 
é"Pli)= D» GX 
k=0 


ce qui ramène notre problème à un problème d’interpolation de Lagrange. 
Remarquons que dès que les c; sont connus, on calcule facilement les a, et 
les b, grâce aux formules 


QAp = Cp FT Cp 


br = Cp — C_»). 
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15.3 Analyse numérique : calcul des intégrales 


15.3.1 Introduction 


Nous voulons calculer numériquement les intégrales définies 


J ds 


où f est une fonction suffisamment régulière pour assurer l’existence des 
dérivées et des bornes supérieures dont nous aurons besoin. 


Les méthodes que nous décrivons dans un premier temps s'appuient sur 
la démarche suivante : 


— On découpe l'intervalle d'intégration en N morceaux de même longueur 


= Ep Li Lis CEy =D 


« — b-a 
OÙ Li+1 — Li — D: 


— Sur chaque morceau on calcule la valeur approchée de l’intégrale 


‘ f(x)dx 


1 


en remplaçant sur l'intervalle [x;,x;,1] la fonction f par une fonction po- 

lynômiale qui interpole f. 
Les diverses méthodes qui entrent dans ce cadre diffèrent par le degré des 
polynômes d’interpolation utilisés et par le choix des points des segments 
[ti, ti1] en lesquels on interpole f. En particulier pour ce dernier point, nous 
verrons Comment choisir au mieux les points d’interpolation. Il faut éviter de 
confondre les points x; du partage en N morceaux du segment initial, avec 
les points des partages des intervalles [x;,x;,1], que nous serons amenés à 
introduire pour appliquer sur ces intervalles des méthodes interpolatoires. 


Comme nous le verrons ces méthodes peuvent être éventuellement complétées 
par une méthode d’accélération de convergence. Enfin nous donnerons un 
abord purement analytique d’une classe de formules de calcul approchée 
d’intégrales. 


102 Chapitre 15 


15.3.2 Mise en œuvre de méthodes interpolatoires 


Dans toute cette subsection pour la clarté des calculs et des méthodes, nous 
travaillerons sur l'intervalle [—1,1] au lieu de [x;,x;,1], puis donnerons les 
formules qu’on obtient de la même façon sur les intervalles [x;,x;,1], et en- 
fin nous sommerons ces diverses formules pour obtenir l’approximation de 
l'intégrale sur [a, b]. Dans chaque cas nous noterons e l’erreur de la méthode 
employée sur [—1,1|, e; l’erreur de la méthode employée sur [x;,x;,1] et enfin 
E l'erreur globale introduite sur le segment [a, b]. 

Signalons enfin que nous supposerons la fonction f à intégrer sur le 
segment [a, b] continument dérivable autant de fois qu’il le faut pour 
pouvoir appliquer les théorèmes permettant l’établissement des majorations 
des erreurs. En particulier on aura à appliquer souvant les théorèmes de 
majoration de la différence d’une fonction et d’un de ses polynômes interpo- 
lateurs. 


Nous noterons 


mP(f)= sup [FPE], 


te[-1,1] 


mP(f)= sup [fP()|, 


ter: ,xi 1] 


MP (f) = sup |fP (|. 


te |[a.b] 


Méthode des rectangles 


La méthode que nous présentons mérite à peine le nom de méthode numérique, 
car elle ne donne pas un calcul efficace et n’est pas employée à cet usage. 
Cependant elle est importante pour la clarté de l’exposé, car elle initialise 
en quelque sorte un processus qui va nous conduire à des méthodes plus effi- 
caces. Grâce à elle nous pourrons présenter et illustrer un certain nombre de 
problèmes de fond. 


La méthode des rectangles, inspirée de la définition même de l’intégrale de 
Riemann consiste à supposer la fonction constante sur tout l’intervalle [—1, 1], 
la valeur de la constante étant celle prise en un point déterminé par la fonc- 
tion. Nous prendrons ici la valeur de la fonction au point —1 (méthode des 
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rectangles à gauche). Remarquons qu’il s’agit bien d’une méthode d’interpo- 
lation de Lagrange en un point par un polynôme de degré zéro. Ainsi 


[ s04- J'rema+e =2f(-1)+e. 


En écrivant 


on obtient 


ou encore 
el < 2m (f). 


Un calcul analogue sur l'intervalle [x;,x;:1] donne 


7 FOdt= (ain af) + es 


avec ? 
Ti+1 — Ti 
lei < mi (9). 
Enfin sur l'intervalle [a, b] nous obtenons la formule globale 


au = = 
[red = rte) +++ flex DI+E 


avec 
@—a) 


lE| < MO(F). 


Méthode du milieu 


La métode des rectangles utilise la valeur de la fonction f en un point 
déterminé de l'intervalle. Il faudrait voir si un bon choix de ce point ne 
conduirait pas à une optimisation de la ” qualité de la méthode”. Encore fau- 
drait il donner un sens à cette expression, ce que nous ferons par la suite. 


104 Chapitre 15 


En attendant exposons la méthode du milieu qui consiste à utiliser comme 
valeur approchée de la fonction sur l'intervalle [1,1], sa valeur au point 
milieu (c’est-à-dire en 0). Ainsi 


[ro - [ ro+e = 2f(0)+e 


En écrivant 


f(u) = #(0) + | fat, 


0 
1 


u f(u)du = 2f(0) + ee Le rot) du. 


Or la méthode d'intégration par partie appliquée aux fonctions 1 et f'(é) 
en dérivant f'(t) et en intégrant 1 (on choisit  — u comme primitive de la 
fonction 1) nous donne 


on obtient 


L | fDat = uf/(0) - Î (tu) FO (Dar. 


1 ( É : (02) ee fl | ( Î - wr%(t) du. 


On en conclut que 


e=|/ EOTE 24(0)] < 1/3 0(9) 


Un calcul analogue sur l'intervalle [x;, x;:1] donne 


EEE Ti + Li 
| FO)dt = (œis1 — 5) T) +e 
avec 5 
Titi — Li 2 
el Em (f). 


Sur l’intervalle [a, b] on obtient par sommation 


roue ir +e+ ras +8 
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où 
ia + Li+1 
Lu 2 
et où , 
b—a) 
El <= ma (f), 


Remarques concernant les deux méthodes précédentes 


Ces deux méthodes sont toutes les deux basées sur une interpolation de degré 
0 (c’est à dire par une constante). Cependant la deuxième a été en quelque 
sorte optimisée en choisissant un bon point d’interpolation. Pourquoi dire que 
cette deuxième méthode est meilleure que la première ? Ceci pour au moins 
deux raisons : d’une part l’erreur E est en 1/N dans la première méthode et 
en 1/N? dans la deuxième, d’autre part la première méthode donne un calcul 
exact pour la classe des polynômes constants alors que la deuxième donne 
un Calcul exact pour une classe plus vaste, celle des polynômes de degré 1 
(ce qui ne veut pas dire qu’au hasard des choix, on ne puisse pas tomber 
sur des fonctions particulières autres, où l’une ou l’autre des formules donne 
une valeur exacte). Nous serons amenés plus tard à préciser ces critères de 
comparaison. 


Méthode des trapèzes 


Passons maintenant à une interpolation de Lagrange de degré 1 (c’est-à-dire 
par une fonction affine) aux bornes de l'intervalle. Ainsi sur l’intervalle [—1, 1] 
on approchera f(u) par 


(D AD Dm D mn) 
Sd 
Alors on aura 


Î f(u)du = f(1) + f(—-1) + e. 


Le calcul de l'erreur e s’effectue en introduisant l’erreur dûe à l’interpolation 
de Lagrange, erreur calculée dans le chapitre portant sur l’interpolation : 


(u + 1)f(1) ce (u — DAT = = 1)(u + 1)!m2 (f). 


fu) - : 
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En intégrant cette inégalité sur [—1,1] on obtient la majoration suivante de 
l'erreur 


el < sm (). 


Des calculs analogues sur l'intervalle [x;, x;:1] nous donnent 


Ju f(u)du = (x) + fait) + ei 


avec ñ 
LP Ve 
a < Mn), 0 (9 


Enfin sur l'intervalle [a, b] on peut écrire 


[rou- Fo) + ON) | 3421) ++ (an) LE 


avec 
@—a) 


ET < 
ET < 12N? 


MA(F). 


Méthode des tangentes 


Cette méthode consiste à approcher la fonction par sa tangente au point mi- 
lieu de l'intervalle. C’est donc une interpolation d’Hermite par un polynôme 
de degré 1. Des considérations géométriques très simples montrent que cette 
méthode se ramène exactement à la méthode du milieu. Aïnsi on dispose 
d’une piste pour expliquer que la méthode du milieu soit meilleure que celle 
des rectangles : en bien choisissant le point d’interpolation de degré 0, tout 
se passe comme si on disposait d’une interpolation de degré 1. 


Méthode de Gauss 


Comme dans le cas de l’interpolation de Lagrange de degré 0, nous allons es- 
sayer dans le cas de l’interpolation de Lagrange de degré 1, de choisir au mieux 
les deux points d’interpolation de manière à améliorer la méthode. Aïnsi pour 
la méthode des trapèzes les point d’interpolation étaient aux bornes de l’in- 
tervalle, pour la méthode de Gauss les points d’interpolation seront des points 
bien choisis de l’intervalle. 
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On se place dans l’espace R.[X] des polynômes de degré < 1 sur R. Soient 
a et 5 tels que —1 < à < 8 < 1. On note 6, et Ôg les formes linéaires sur 
R:[X1] définies par 0,(P) = P(u). Si on note 


XN — 
PA) = TS 
et > 
Ps(X) = — 


alors d(P;) = Ou (0 si u À v, 1 si u = v). On en conclut que (ô4, 05) est 
une base du dual de R,[X]. Donc la forme linéaire 7 su R[X] définie par 


1 
I(P) =) f(u)du 
ci 
se décompose sous la forme 
PASSE À808. 


Cette formule appliquée aux polynômes P, et P; donne les valeurs explicites 
des constantes À, et Àg, ce qui permet d’écrire 


28 2a 
Eee 


Ceci veut dire que pour tout polynôme P € R;[X] on a 


T'= 


ôg. 


| _—.  — 
Le P(u)du = re re 


Pour une fonction f on va donc avec cette méthode écrire 


[fou = ro) + HO +e 


et l’erreur e est nulle pour les polynômes de degré < 1. Peut-on choisir 
les points & et B pour que la formule reste exacte (e = 0) pour un degré 
supérieur ? La réponse est donnée par : 
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Proposition 15.3.1 1! existe un couple et un seul de points a et B tels que 
pour tout polynôme P de degré < 3 on ait 


[roue Pt + pie 


ces points sont 
—1 8 1 
a = — Re 
V3 V3 


De plus, quels que soient les points a et 5 il existe un polynôme de degré 4 
qui ne vérifie pas la formule précédente. 


Preuve : le calcul se fait à partir de la formule exacte pour les polynômes 
de degré < 1 en imposant que cette formule reste exacte pour le polynôme 
X? et le polynôme X* (ce qui est nécessaire et suffisant pour que la formule 
reste exacte pour les polynômes de degré < 3). Un calcul simple nous donne 
les points à et 6 attendus aïnsi que les coefficients explicites de la formule. 


Proposition 15.3.2 Pour tout polynôme de degré < 3 on à 


: 1 1 
P(u)du = P(—-—>) + P( 
fe 7 
La méthode de Gauss consiste donc à utiliser cette formule pour approcher 
l'intégrale : 


1 
1 1 

Fudu = f(-) + FS 
É VAL V3 
La majoration de l’erreur e se fait maintenant en utilisant un polynôme 


de degré < 3 qui interpole f. Plus précisément on considère le polynôme 
d’'Hermite H de degré < 3 tel que 


)+e. 


1 l 1 1 
HE) AC AIS) 
1 L LL gl ue Ta … re 1 
CE ICT) Hs) = 5) 
Dans ces conditions 
À 1 1 
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e=|/ EURO TE ROEON 


On sait alors (cf. chapitre sur l’interpolation) que sur [—1,1] 


donc 


1 1 
= Lea ee (4) 
fu) — Hu) < 5-3) m(p) 
d’où on tire par intégration 
(4) 
a < 90). 
135 


Sur chaque intervalle [x;,x;:1] on obtient de manière analogue 


2 2 


à F(é)dt = Ti+1 — Li E | — + Ti je Ti+1 — “) +f (2 + TL 


2 23 


avec 


4 
le;| < Li 2) m 2) 


Par sommation on obtient sur [a, b] 


b NET pi: 
| f(u)du = > f(u)du + E 
a i=0 * 


avec 
a) MO(F) 
NT 4320 


LE < 
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Remarque : En fait bien qu’au départ on ait travaillé sur une interpolation 
de Lagrange de degré 1, tout se passe ensuite de la même façon que pour 
une interpolation d'Hermite de degré 3. C’est ce qui fait la performance de 


la méthode. 


Méthode de Simpson 


Nous partons maintenant de l’interpolation de Lagrange de degré 2 (par une 
fonction parabolique) aux bornes de l’intervalle et en un point 7 tel que 
—1 < y <1. Nous employons une démarche analogue à celle du paragraphe 
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précédent. Notons 6_:1,0.,01 les formes linéaires sur l’espace R2[X] des po- 
lynômes de degré < 2 définies par 


ôa(P) = P(u) 
Les polynômes ee (X —1(X — +) 
2(1+79) 
PE RE 
AO = STE 


vérifient 0,(P,) = ôuy. On en déduit que (8_:,9,,91) est une base du dual de 
R2[X] et donc que la forme linéaire Z définie par 


se décompose sur cette base 
I — À_10 1 + À:07 + M1. 
Ceci veut dire que pour tout polynôme P de degré <2 on a 
: 
| P(bt = À 1P(-1) + À, P(9) + MP). 
= 


Par un bon choix du point 7 on va voir que cette formule persiste pour les 
polynômes de degré < 3. On laisse au lecteur d’établir la proposition suivante 


Proposition 15.3.3 Il existe un point y et un seul tel que pour tout po- 
lynôme de degré < 3 on ait 


. P(t)dt = À_1P(—1) + À, P(7) + ÀP(1). 


1 


Cette valeur de 7 est 0. De plus il existe un polynôme de degré À pour lequel 
cette formule n’a pas lieu. 


Ce point y = 0 étant choisi on obtient la formule dite des 3 niveaux. 
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Proposition 15.3.4 Pour tout polynôme de degré < 3 


| | PA = IP(-1) + P(1) + 4P(0)]. 


1 


La méthode de Simpson consiste donc à utiliser cette formule pour ap- 
procher l'intégrale : 


L 
1 
Î ftodu = 1-0 + F0 +401 + € 
1 
Soit H le polynôme de degré < 3 tel que 


H(-1)=f(-1) H(0)= (0) H()=f() H(0)= f(0). 


On a alors 


jé [ be [tou 


On sait (cf. chapitre sur l’interpolation) que sur [-1,1] 


SOON 


fu) — Hu] < que = D(u + Dim() 


ce qui donne par intégration 


Sur les intervalles [x;,x;,1] on obtient 


avec 


Ti+1 — Li Ù ms” (f) 
lei < | ————— 
2 90 


Sur l'intervalle [a, b] tout entier on peut écrire 


[so = _. te + fe) +29 194493 es 


i=1 i=0 
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avec 
Ga) MO(F) 
NT 2880 


L'erreur est là aussi comme pour la méthode de Gauss en 1/N#, mais avec des 
coefficients un peu moins bons. Cependant la formule ne fait pas intervenir 
comme dans celle de Gauss des nombres compliqués” (43) ce qui explique 
qu’elle ait été plus prisée avant l’apparition des ordinateurs. 


|E| < 


Peut on là aussi optimiser mieux la méthode en choisissant de manière astu- 
cieuse les trois points d’interpolation ? Ceci est en partie l’objet de l’étude 
plus générale du paragraphe suivant. 


Méthodes interpolatoires et polynômes orthogonaux 


Soient P5, P1,---,P, les n + 1 premiers polynômes de Legendre. Rappelons 
que les polynômes de Legendre forment l’unique suite de polynômes telle que 
a) deg(P;) = i 
b) fa Pi(u)P;(u)du=0 (4j) 


c) P() = 
Appelons a,---,a, les racines de P, : elles sont distinctes et appartiennent 
à l'intervalle [-1,1]. Les formes linéaires (ô4.,---,0,,) définies par ô4,(P) = 


P(a;) forment une base du dual de l’espace R,;_1[X] des polynômes de degré 
< n—1 (espace de dimension n). On en conclut que la forme linéaire J définie 
sur R;_1[X] par 


se décompose sur cette base sous la forme 
I — Aa 0e, + Ut + AO 


ce qui veut dire que pour tout polynôme P de degré <n—1ona 
1 
Î Ptudu = x Pa) +++: + 2, P(an) 
4 


Le choïx des racines du polynôme de Legendre P, pour les points a; est un 
très bon choix dans le mesure où 
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Proposition 15.3.5 Pour tout polynôme de degré < 2n — 1 on a 


Î P(u)du = a, P(a) +++ 4, P(an). 


1 


Preuve : Il suffit de montrer que la formule est vraie pour les polynômes 
d’une base de l’espace R2,_1[X] des polynômes de degré < 2n — 1. La famille 
(Po, - +, Pn_1; PoPns* **, Pn_1Pn) est une base de R2,_1[X] (tous les degrés 
sont représentés une fois et une seule). Il est clair que Po,---, P,_1 vérifient 
la formule puisque ces polynômes sont de degré < n — 1. Quand aux P;P, 
(0 <i<n—1l)ilest clair en utilisant l’orthogonalité de la suite des polynômes 
de Legendre que [', P(u)P,(u)du — 0 et que sachant que les a; sont les 
racines de Ph, 6(P;Pn) = Pj(ai)Pa(ai) = 0. Ils vérifient donc eux ausi la 
formule. 


Toute autre famille de points est moins bonne que la famille (a;). 


Proposition 15.3.6 Soient a:1,--:,an des points distincts de l'intervalle 
[1,1] dont l’un au moins n’est pas racine du polynôme de Legendre de degré 
n. Alors il existe un polynôme Q de degré < 2n — 1 tel que 


de Q{u)du £ An Q() +--- +20, Q(an). 


Preuve : Posons Q,(X) = (X — æ&)---(X — a,).Par hypothèse, Q, n’est 
pas proportionnel au polynôme de Legendre P,, il n’est donc pas orthogonal 
au sous espace des polynômes de degré < n — 1. Il existe un polynôme R 
de ce sous espace tel que fE: R(u)Qh(u)du Æ 0. Il suffit alors de prendre 
Q = RQ. 


Cette propriété optimale des racines des polynômes de Legendre peut être 
utilisée pour avoir une meilleure majoration de l’erreur quand on remplace 
l'intégrale à calculer par l’intégrale du polynôme d’interpolation de Lagrange 
de degré < n—1 en n points de l'intervalle. C’est ainsi que nous avons procédé 
pour la méthode de Gauss, avec n = 2. Les point —1/V3 et 1/13 sont les 


racines du polynôme de Legendre de degré 2 (P(X) = 1), 


Dans le cas de l’interpolation de degré 2 (avec les notations précédentes 
n = 3), nous pouvons obtenir mieux que par la méthode de Simpson et par la 
méthode de Gauss. La méthode de Simpson en effet est une méthode utilisant 
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l’interpolation de degré 2, partiellement optimisée par le choix du point milieu 
de l'intervalle, mais qui n’est pas complètement optimisée puisque les deux 
autres points d’interpolation sont les bords de l’intervalle et non pas des 
racines du polynôme de Legendre de degré 3. Le polynôme de Legendre de 
degré 3 est 


5X?—3X 
P3(X) = DE 
ce qui donne pour racines 
v3 " V3 
UM=-——= à A3 = —=. 
UV V5 


Un calcul simple donne alors pour tout polynôme P de degré < 5 


Î P(u)du = . se (2) + 8P(0) +5P (2) : 


Ainsi pour toute fonction f de classe CŸ on a 


Î f(u)du = . 5P (-2) + 8P(0) +5P (à 


+e 


V5 V5 


où |e| se majore en utilisant l’inégalité 


fu) = P(u)| < Sul — 3/5) 0 (9) 


ce qui par intégration fournit 


1 
€ m6)(f), 
el < sm" ) 


Sur chaque intervalle [x;,x;,1] on à une erreur 


Ti4+1 — Li | m9 (f) 
EC < 3 
2 15750 


et sur [a, b] une erreur 


a) MO(F) 
NS 15750 x 128 


LE] < 
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15.3.3 Présentation générale des quadratures élémentaires 
- Ordre d’une méthode 


Toutes les méthodes vues jusqu’à présent entrent dans le cadre général sui- 
vant (méthodes de quadrature élémentaires). On a sur l'intervalle [x;, x;,1] 
une formule du type 


Titi li 
/ ftu)du = Gun — 25) um) + es, 


Ÿ= 
où 
li 
mis Elrotinl] et ) 5, = 1. 
j—0 


On essaie d'adapter les paramètres pour que la méthode soit la meilleure 
possible. 


Définition 15.3.1 On dit qu’une méthode de quadrature est d’ordre N si 
la formule approchée est exacte pour les polynômes de degré < N et inexacte 
pour au moins un polynôme de degré N +1. 


: L , . : 
Remarquons que puisque }};_çur,; = 1, une méthode de quadrature élémentaire 
est toujours au moins d’ordre 0. 

Le lecteur est invité à préciser les paramètres qui correspondent aux diverses 
méthodes décrites précédemment. 


15.3.4 Accélération de convergence - Formule de Ri- 
chardson - Méthode de Romberg 


La méthode que nous allons décrire maintenant part d’une méthode d’accélération 
de convergence sur les suites. 


Formule de Richardson 


Soit u une suite de nombres réels convergeant vers un nombre réel a. On 
suppose que 


Un — a = ÀK" +O(k") 
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où k et k’ sont des réels tels que 0 < |k'] < |k| < 1 et où À est un réel non 
nul. 
Alors 

Una — a = AK + O(k") 


kun — ka = Xk"*1 + O(k") 


et par suite 
Un+1 — kun 


1—k 
On construit de cette manière une suite 


— a = O(K"). 


de Un+1 — kun 
" 1—k 
qui converge vers à plus vite que la suite initiale u. 
Plus généralement on peut supposer que 
où 
0 <kpul <lkl <<] <1 


et où ÀA1,---, À, sont des réels non nuls. Pour toute suite w, notons R£w la 
suite définie par 
Un+1 — EWn 


RU —= 1— E 
Alors k 
Un+1 — Ain 
Run = — — 
ki 4 1—k 
Puisque 
Uni Q= À kT TT + Ok + eee + ARE + O(KP. 1) 
et que 
KiUn = ka —= MR + Ak1ks His Àkik, + O(kr11) 
on obtient 


Un+1 — kiUn = (a = k1a) —= (AkS ++ À,k7)( = k1) + O(kr11) 
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où 
(ko — ki) 


1 — ... —= —————— 
co LR re PER 


On en déduit que 
Rain — a = À ++ AK) + O(krii). 
Par conséquent on peut itérer le procédé et on obtient alors 
Ry,Rpo 1" Rriln — à = O(kr). 


Remarquel : Reprenons la formule de Richardson à l’ordre 1. Il se peut 
que k ne soit pas connu. On pose alors 


Un+1 — Unün 
Ve = 


ENT 
où 
Un+1 — Un 
Un = ——— 
Un — Un-1 


Dans ces conditions 


Un —k=O((K/E)), 


donc 

Dao =O(ET) 
Remarque 2 : On peut donner une version continue de la formule de 
Richardson. 


Soit f(y) telle que lim,_,0 f(y) = &. On suppose que 
fu) = à + ay ++ ay + O(yTt). 


Soit alors 0 < r < 1 et yo > 0. 
Formons la suite 


Arm,0 — f(r""v), 
et remarquons que 


lim Aym,0 — 0, 
mMm—00 


Am, = O0 + @ÿor® +. + ar?) + O ((r7#7)7). 
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Donc en posant ki =r,---,k, = 7? on se ramène à la formule de Richardson 
précédente. Plus précisément, pour 0 < n < p — 1 on définit 


+1 
ns ut An 


Apart = 1 —rv+#1 


Dans ces conditions 
Amp — %0 = O (éd ie) : 


Méthode de Romberg 


Nous allons montrer comment marche cette méthode . 
Découpons l'intervalle [a, b] par un partage équidistant en N = 2” morceaux 
et posons À = (b — a)/2". Appelons 


27 —1 


In = 1/2h{f(a) + (6) +2 3° F(x:)] 


la valeur approchée de l’intégrale donnée par la méthode des trapèzes. 
La formule d’Euler-Maclaurin nous donne 


[l f{u)du — In = a(1/4)" + O ((1/8)"). 


On peut donc appliquer la méthode de Richardson à l’ordre 1 avec k = 1/4 et 
k' = 1/8. On pourra remarquer que dans ce cas on obtient la formule donnée 
par la méthode de Simpson sur 2°! intervalles. Il faut appliquer la méthode 
de Richardson à un ordre > 3 pour trouver une formule qui n’est pas donnée 
par les considérations des paragraphes précédents. 


Bibliographie 


[1] J.P. Demailly Analyse numérique et équations différentielles 
Presses Universitaires de Grenoble 1991 


[2] I. Glazman, Y. Liubitch Analyse Linéaire dans les Espaces de 
Dimensions Finies Mir-Moscou 1972 

[3] M.W. Hirsh, S. Smale Differential Equations, Dynamical Sys- 
tems, and Linear Algebra Academic Press 1974 

[4] J.L. Ovaert, J.L. Verley Algèbre Vol.1 Cedic Nathan (Coll. Epis- 
temon) 1981 

[5] R. Rolland Groupes finis commutatifs et transformations de 


Fourier discrètes Cours de L’ES2I 1992 (disponible par ftp sur 
iml.univ-mrs.fr) 


119 


Index 


accélération, 117 

algèbre de Boole, 9 
algorithme de Crout, 63 
application linéaire, 15, 18, 55 


base, 20, 23, 55 
base duale, 34 

base incomplète, 24 
bidual, 34 

bloc de Jordan, 77 
Boole (algèbre), 9 


caractéristique (polynôme), 67, 75 
caractère, 80 


Cayley Hamilton (théorème de), 75 


circulante (matrice), 43 
codimension, 28 
cofacteur, 57 
combinaison linéaire, 19 
convolution, 81 

Cramer (système de), 62 
Crout (algorithme de), 63 


décomposition spectrale, 68 
dépendant, 20 

déterminant, 50 

déterminant de Vandermonde, 53 
diagonal, 74 

différences divisées, 85, 87 
dimension, 24 

dimension finie, 23 


120 


directe (somme), 27, 30 
droite, 25 

dual, 33 

duale (base), 34 


échange, 23 
élimination, 63 
endomorphisme, 15 
engendré, 12 
espace engendré, 12 
espace quotient, 56 
espace vectoriel, 7 


factorisation des applications, 56 


famille génératrice, 20, 23 
famille liée, 20 

famille libre, 20, 23 

filtre stationnaire, 81 

forme linéaire, 18, 33 

forme multilinéaire, 49 
formule des trois niveaux, 35 
Fourier (transformée de), 91 


générateur, 20, 23 
Gauss (méthode de), 108 
Gauss (pivot de), 63 


Hermite (interpolation d’), 96 
homomorphisme, 15, 18 
hyperplan, 25, 34 


image, 16 


incomplète (base), 24 
indépendant, 20 

index, 69 

interpolation, 82 

interpolation d’Hermite, 96 
interpolation de Lagrange, 83, 91 
inverse (matrice), 57 


Jordan (bloc de), 77 


Lagrange (interpolation de), 83, 91 
Lagrange (polynôme de), 84 
Legendre (polynômes de), 114 
liée, 20 

libre, 20, 23 

linéaire (application), 15, 18 
linéaire (combinaison), 19 
linéaire (forme), 18, 33 
linéaire (système), 59 
linéairement dépendant, 20 
linéairement indépendant, 20 


méthode de Gauss, 108 
méthode de Richardson, 117 
méthode de Romberg, 117 
méthode de Simpson, 111 
méthode des rectangles, 104 
méthode des tangentes, 108 
méthode des trapèzes, 107 
méthode du milieu, 105 
matrice, 40, 45 

matrice circulante, 43 
matrice de passage, 56 
matrice inverse, 57 

matrice transposée, 46 
milieu (méthode du), 105 
minimal (polynôme), 69 
multilinéaire (forme), 49 


121 


Newton (polynômes de), 84 
noyau, 16 


opérateur, 15 
ordre d’une méthode, 116 
orthogonaux (polynômes), 114 


passage (matrice de), 56 

pivot de Gauss, 63 

polynôme caractéristique, 67, 75 
polynôme minimal, 69 
polynômes de Lagrange, 37, 84 
polynômes de Legendre, 114 
polynômes de Newton, 84 
polynômes orthogonaux, 114 
produit de matrices, 45 
projecteurs spectraux, 74, 76 
projection, 28, 31 

propre (sous espace), 68, 70 
propre (valeur), 67, 70 

propre (vecteur), 67 


quotient (espace), 56 


rang, 24, 55 

rectangles (méthode des), 104 
Richardson (méthode de), 117 
Romberg (méthode de), 117 


signal, 79 

Simpson (méthode de), 111 
somme, 27, 30 

somme de matrices, 45 
somme directe, 27, 30 

sous espace engendré, 12 
sous espace propre, 68, 70 
sous espace spectral, 70 
sous espace vectoriel, 11 
spectral (projecteur), 74, 76 


122 


spectral (sous espace), 70 
spectre, 70 

splines, 98 

stationnaire (filtre), 81 
supplémentaire, 28, 56 

système de Cramer, 62 

système de Vandermonde, 84, 88 
système linéaire, 59 


tangentes (méthode des), 108 
théorème de division, 86 
transformée de Fourier, 91 
transposé, 38, 46 

trapèzes (méthode des), 107 


valeur propre, 67, 70 
Vandermonde (déterminant de), 53 
Vandermonde (système de), 84, 88 
vecteur propre, 67 

vectoriel (espace), 7 

vectoriel (sous espace), 11 


